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Preface 


This is a book on empirical number theory concentrating on the analysis of 
number sequences. Its focus is on a small part of a very large number of 
integer sequences defined by Florentin Smarandache. The author has, 
however, when appropriate included some other of his research results which 
organically belongs to this area. The content is organized into chapters 
according to the main considerations to make when programming the 
analysis. They are not mutually exclusive. 


As in my previous book Surfing on the ocean of numbers - a few 
Smarandache Notions and Similar Topics an attempt has been made to 
present results so that they are easy to understand and not to burdensome to 
read. There are many tables, some of which may be used for reference, but in 
general they are there to show the overall results obtained at a glance. In most 
cases it is the way in which sequences behave, not the individual figures, 
which is of importance. Some graphs have been included to illustrate 
important findings. 


Some of the results in this book were presented by the author at the First 
International Conference on Smarandache Notions in Number Theory, 
August 21-23, 1997, Craiova, Romania. 


References have been given after some of the chapters. However, constant use 
has been made of the following Smarandache source materials: 

Only Problems, Not Solutions!, Florentin Smarandache. 

Some Notions and Questions in Number Theory, C. Dumitrescu, V. Seleacu. 


Illustrations, graphics, layout and final editing up to camera ready form has 
been done by the author. Most tables have been created by direct transfer from 
computer files established at the time of computing to the manuscript so as to 
avoid typing errors. The results often involve very large numbers which are 
difficult to accommodate in pre-designed formats. It has therefore often been 
necessary to use several lines to represent a number in tables as well as text. 
ΑΗ calulations have been carried out using UBASIC, ver. 8.87. 


I am grateful о Dr. R. Muller who has given all possible help and 
encouragement during the work on this book. He and his colleagues at the 
American Research Press have at all times facilitated the work through rapid 
e-mail communications. 


Finally, it is thanks to the patience and understanding of my wife Anne-Marie 
that this book has come about. But the winter months have passed, the book is 
finished and I am no longer going to be lost among manuscript pages. 
Summer is around the corner and the Swedish nature is waiting for us. 


Paris, April 1998 
Henry Ibstedt 
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Chapter I 


Partition Sequences 
L1 Introduction 


We will study the partition of the sequence of digits 
of positive integers into two or more groups of digits whose corresponding 
integers may or may not possess a certain property. 
Example: 
17|7|101 
Is a partition of the 6-digit integer 177101 into three groups of digits (parts) 
each of which represents a prime, whereas 
1|77|10|1 
is a partition of the same integer into four parts none of which is а prime. 


Obviously there are many ways in which and integer can be partitioned into 
groups of digits. The kind of partition described above may also be referred to 
as a de-concatenation in order to distinguish it from the classical concept of 
partition of an integer as first developed by Euler and later illustrated by 
Ferrer [1]. However, the kind of partition we are dealing with will be clear 
from the context. We will now use the classical partition of an integer as a 
means to help us formulate a strategy for partitioning sequences of integers in 
the sense defined above. 

The partitions of 4 are: 1+1+1+1=2+1+1=2+2=3+1=4 with the corresponding 
Ferrer diagram: 


OOOO 0O UU OO O 
لا‎ OO LJ 


ev JE] 


1 2 3 4 


Diagram 1.1116 Ferrer diagram for 4. 


The last partition corresponds to 4 itself which is of no interest to us. 
Furthermore, we will link the vertical representations to form the pattern 
shown in diagram 2: 

ΠΠ COO LCA Coo 


1 2 3 4 


Diagram 2. The Ferrer diagram for 4 rearranged and with partition numero 5 
excluded. 


In our applications we will also need to consider the order in which the 
partition elements occur. Diagram 2 shows the seven different ordered 
partitions possible for 4 and how they form a partition pattern for the 
partition of the 4-digit integer 9164 into groups of integers. 


(9121614) [si (6/41 [136 Y«] Sie) [91] 6«] [51e [4] [J τες] 
1 2 3 4 5 6 7 


Diagram 3. The partition pattern for а 4-digit integer. 


From this we see that the integer 9164 can be partitioned into perfect squares 
in exactly two ways, 9, 16, 4 (partition numero 3) and 9, 1, 64 (partition 
numero 4). This is how the mapping of and integer onto a partition pattern 
helps us study the properties of each partition element. 


Only partitions 2 and 4 in diagram 2 contain partition elements of unequal 
sizes. These give rise to an increase of the number partition patterns through 
arrangement of the partition elements. It is easily understood that in a 
partition of 8 we have 2 parts of the type لا‎ and 4 parts of the type O 
then the number of arrangements will be 2415. The integer 8 has 21 
partitions of the type shown in diagram 2 giving rise to 127 ordered 
partitions. This could be calculated by considering each arrangement as we 
did above. However, we will soon arrive at this in a different way. 





The above considerations will be useful in order to express the partitions of 
integers in explicit form. When we consider the partition of a very large 
number of integers to find out how many of them possess a certain property 
we will be able to use a different approach. 


Let's represent an n-digit integer t in the form t=a,+a,-10+a,-107+ ... а,.10°!. 
We introduce the following definitions and computer analysis related 
concepts in UBASIC. 


O(t}=1 ift possesses a certain property otherwise 5(t)=0. 
\ symbolizes integer division, example 23678\100=23 


res symbolizes the remainder of the last performed integer 
division, in the example above res=678. 


We will now take a look at the question of leading zeros. Is 3 | 07 a partition 
of 307 into two primes? Is 8 | 027 a partition of 8027 into two cubes? The 
author prefers to have a unique representation of an integer, 7 is a 1-digit 
integer written 7 not 07 and 37-27 is written 27 - not 027. The computer, 
however, interprets 07 as 7 and 027 as 27. To avoid integers with leading 
zeros the function е(),г) is defined as follows. 


g(j,r)=0 if r<10"' otherwise e(j,r)=1, where г and j are context 
depending variables. 


We will now use these functions to formulate an algorithm for calculation of 
the number of n-digit integers which can be partitioned so that the partition 
elements have a given property (being primes, squares, cubes, etc) 

Lets denote the number of partitions of t by p(t). 

l-digit integers: 1<t<9. Obviously there is nothing to partition, p(t)-0. 


2-digit integers: 10<t<99. p(t) can only assume the values 0 and 1 and is 
given by 
p(ty=S(t10)-6(res) 


3-digit integers: 100<t<999. In the algorithm below use is made of the fact 
that we know p(t) for t«100. 


p(t)-5(t100)-e(1,res)-((res)* p(t100))5(00)-(5(res)* p(t10)), 


where we know that p(110)-0. It has been included in order to present a 


n-digit integers: 10*' «t«10"-1. 


n-1 


P(t) = > Ot \10"*) (п - k —1,resXó(res) + p(res)) (1) 


ΕΞ] 


We сап now use (1) to calculate the number of partition patterns for n-digit 
integers. To do this we put S(t\n™“}=1 (the property of being a partition 
element), e{n-k-1,res}=1 (a partition element may begin with а zero) and 
p(res)-p(k). This results in 


п-1 
p(n) = 2 (1+ p(k)) , where p(1)-0 
k=! 
Evaluation of the recursion formula gives p(n)-2"'-1 resulting in the 
following table of partition patterns. 
Toble 1. Partition pattems, N = number of pattems. 


n-digit integer } 2 3 4 5 6 7 8 
М 0 L 3 7 15 3l 63 127 


The algorithm has been implemented іп UBASIC оп a Pentium 100 Mhz 
portable computer to examine various partition sequences. 


1.2 Тһе Smarandache Prime-Partial Digital Sequence 


Definition: The prime-partial-digital sequence is the sequence of prime 
numbers with the property that each member of the sequence can be 
partitioned into groups of digits such that each group is a prime. 


Table 2 shows the first one hundred terms of this sequence. 


We will see in Chapter Π.4 that the Smarandache Prime-Digital Sub- 
Sequence is infinite. The Prime-Partial Digital Sequence is in fact a subset of 
the Prime-Partial Digital Sub-Sequence and is therefore a fortiori infinite 
which settles affirmatively a conjecture by Sylvester Smith. [2]. We will make 
Observations on the number of ways in which prime integers can be 
partitioned into primes. How close is the maximum number of ways this сап 
be done to the theoretical maximum given in table 1 and what do these 
partitions look like? 


Tabie 2. The first one hundred members of the prime-partiat-digit sequence. 


23 37 53 73 113 137 173 193 1987 211 
223 227 229 233 241 257 271 277 283 293 
311 313 317 331 337 347 353 359 367 3 
379 383 389 397 433 523 541 547 557 571 
577 593 613 617 673 677 719 727 733 1743 
757 761 773 797 977 1013 1033 1093 1097 1117 

1123 1129 1153 1171 1277 1319 1327 1361 1367 1373 
1493 1637 1723 1733 1741 1747 1753 1759 1777 1783 
1789 1913 1931 1933 1973 1979 1993 1997 2113 2131 
2137 2179 2213 2237 2239 2243 2251 2267 2269 2273 


The property that we will examine with the function (0) is whether t is а 
prime ог not. In UBASIC this is done іп the following program. 


e(t): 10 .D(t) * D(t) is equivalent to δ(8) 
20 if t-0 then 2-0 :goto 40 
30 if nxtprm(t-1)-t then z-1 else z-0 
40 return(z) “7 is the value of D(t) 
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The function e£(j,r) takes the following form in UBASIC. 


e(j,r): 10 LG% r) ўт) is equivalent to 
e(j,r) 

20 if r<10°G%-1) then y=0 else y-1 

30 return(y) «y is the value af e(j,r) 


These two functions were made work together with the algorithm (1) in a 
UBASIC program to produce the results shown in tables 3 and 4. 


Table 3. The ratio т/а where m is the number of n-digit primes which can be 
= Ῥσποηπθα into primes ond q is the number of n-digit primes. 


^, 





"1-2 |^ n=3 =4 =S n=6 و‎ nz8 
m. 4 2 51 383 3319 27111 234229 2040170 
а 2X ^ 143 1061 8363 68906 586081 5096876 


m/q 0.1908 0.3566 0.3610 0.3969 0.3934 0.3997 0.4003 








Diagram 4. The ratio т/а for n-digit integers (n—2,3, ...8). 
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Table 4. The number of n-digit primes that can be partitioned into primes in p 
different ways. (The tail for n=8 is interposed.) 








0 17 92 678 504 41795 351852 3056706 
| 4 40 251 1986 14789 118379 972095 
2 10 83 769 6306 54102 436009 
3 1 40 309 2811 25892 228428 
4 5 135 1445 14205 139755 
5 4 67 758 7676 5 
6 33 424 5012 55793 
7 15 250 2922 35326 
8 3 131 1877 24401 
9 1 78 1258 17131 

10 0 46 873 12424 

1 0 22 550 8724 

12 1 27 481 6671 

3 : p n-8 ; 12 277 4841 

14 : 38 16 i 6 212 3573 

15 : 3۶ 17 ἱ 3 147 278] 

16 : 40 8 i 0 115 2052 

7 41 ІЗ! 0 62 1574 

18 42 5i 3 59 1244 

19 43 4i 33 992 

20 44 8: 34 811 

21 45 74: 18 617 

22 46 3i 12 500 

23 47 9: 6 417 

24 48 ό | 10 311 

25 49 2 4 264 

26 50 2i 4 214 

27 5] 4 | 2 154 

28 52 / | 3 150 

29 53 0: 0 94 

30 54 4: | 83 

31 55 0: 2 85 

32 56 2 0 65 

33 57 ο i 0 5] 

34 58 0 i 0 42 

51 0 0 وی 35 

36 60 ΔΕ 0 30 

37 1 32 
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Table 5. All possibie partitions of 23733737 into primes. 
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Table 5. Continued. 


42 2 37 3 373 7 
43 2 3 73 373 7 
44 2 373 373 7 

45 2 37337 3 7 

46 2 337 37 37 

47 23 733 73 7 

48 23 733 7 37 

49 23 7 337 37 

50 2 37 337 37 

51 23 73 373 7 

52 2 3733 73 7 

53 2 3733 7 37 

54 23 7 3373 7 

55 2 37 3373 7 

56 2 37337 37 

57 23 3 7 

58 23 7 33737 

59 2 37 33737 

60 2373373 7 





With a slight exception for 6-digit integers the ratio between primes that can 
be partitioned into primes and the number of primes is increasing and 
indicates strongly that the prime-partial digital sequences is infinite. 


To list the actual partitions would be impossible but for each n-digit integer 
sequence of primes there are a few primes that are record holders in the 
number of ways in which they can be partitioned into primes (p ways). Here 
is a list of those primes. 


n=3 p=3 373 

n=4 p=5 3137, 3373, 3733 and 3797 
n=5 p=12 37337 

n=6 p=18 237373, 537373 and 733373 
n=7 p=37 2373373 

n=8 p=60 23733737 


We see that of the prime digits 2,3,5 and 7 the primes 3 and 7 play a 
dominant role in the composition of those primes that lend themselves to the 
maximum number of partitions into primes. All partitions can be displayed in 
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explicit form by executing a UBASIC program based оп the partition patterns 
which were subject to a detailed discussion in the introduction. Table 5 
shows the partitions for 23733737. 


L3 The Smarandache Square-Partial Digital Sequence 


Definition: The square-partial digital sequence is the sequence of perfect 
squares which can be partitioned into two or more groups of digits which are 
also perfect squares. 


The first one hundred terms of this sequence is shown in table 6. 
Tobie ὁ. The first one hundred terms f the square-parfial digital sequence. 


49 100 144 169 36] 400 441 900 125 139 
1444 1600 168) 19% 2500 329 3600 425 4900 6400 
8100 9025 9409 10000 10404 11025 11449 11664 12100 12544 
14161 14400 14641 15625 16641 16900 19044 19600 22500 25600 
28900 32400 36100 36481 40000 40401 41616 42025 43581 44100 
44944 48400 49729 52900 57600 62500 64009 67600 72900 78400 
81225 84100 90000 93025 93636 96100 99225 102400 105625 108900 

115600 116964 117649 119025 121104 122500 129600 136161 136900 140625 
143641 144400 152100 157609 160000 161604 164025 166464 168100 170549 
176400 184900 193600 194481 202500 211600 220900 225625 230400 237169 


This sequence is infinite. If all the infinitely many squares of the form 5.103, 
where s is a perfect square, are removed the sequence is still infinite. Charles 
Ashbacher [3] has proved this in the case where leading zeros are accepted. 
In order to continue to use unique representation of integers the proof has to 
be changed. 


Proof: There are infinitely many squares of the form (10334 2). In expanded 


form these can be written (105/7.107**9.«(2.1057.10 +27. From this we see 
that (10%*'+2)* can be partitioned into the three squares: 10%, 4-10% and 4. 
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The methods and programs developed in the introduction to this chapter have 
been used in the analysis of this sequence with the only change that we study 
a different partition property defined though the function δίῃ which now 
takes the form: 


10 .D(t) 
20 if (isqrt(t))^2-t then 2-1 else 2-0 
30 return(z) 


Use has been made of the Ubasic function isqrt(x) which gives the integer 
part of the square root of x. 


The analysis was carried out for 2, 3, ... 8 digit squares. The results are 
presented in the form used in the previous section. 


Table 7. The ratio m/q. mis the number of squares which can be partitioned into 
squares and q is the number of n-digit squares. 


n=2 n=3 n=4 n=5 n=6 п-7 n=8 


m 1 7 15 44 114 316 883 
q 6 15 68 217 683 2163 6837 


m/q 0.1667 0.4667 0.2206 0.2028 0.1669 0.1461 0.1292 





Diagram 5. The ratio m/q for n-digit integers (n = 1,2 ... 8). 


For squares with more than three digits the graph shows an interesting 

decrease as if it were approaching an horizontal asymptote. Is this the case? 

Doe ο πο 
be partitioned into squares in a maximum number of ways (p ways): 


n=3 p-1 100, 144, 169, 361, 400, 441, 900 
n=4 р-3 4900 

-5 p=5 11449, 14400, 16900, 36100, 44100 
n-6 . -5 144400, 490000 

n=7 p= 1144900 

n=8 p-8 14440000, 36144144, 49491225 


This list is dominated by squares ending on an even number of zeros. The 
digit 4 occurs most frequently while the digits 7 and 8 are lacking. This 
section 15 terminated by a table of maximum partitions for n=8. 


Table 8a. All possible partitions of 49491225 into squares. 


1 4 9 4 9 1 225 
2 4 9 4 9 1225 
3 49 4 9 1 225 
4 4 9 49 1 225 
5 49 49 1 225 
6 49 4 9 1225 
7 4 9 49 1225 
8 49 49 1225 
Table 8b. Al itions of 36144144 into s es. 
1 36 1 4 4 1 4 4 
2 361 4 4 1 4 4 
3 36 144 1 4 4 
4 36 1 441 4 4 
5 36 1 4 4 144 
6 361 441 4 4 
7 361 4 4 144 
8 36 144 144 
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Table 8c. All possible partitions of 144440000 into squares. 


1 1 4 4 4 0 0 0 0 
2 144 4 0 о 0 0 

3 1 4 4 400 0 0 

4 1444 0 0 0 0 

5 144 400 0 0 

6 1 4 4 40000 

7 144400 0 0 

8 144 40000 





L4 The Smarandache Cube-Partial Digital Sequence 


Definition: The cube-partial digital sequence is the sequence of cubes which 
can be partitioned into two or more groups of digits which are also cubes. 


It is obvious that all integers of the form n°.10™ where n is an integer belong 
to the sequence, but if we let n be an integer whose last digit #0, then there 
are only two cubes m-n'«10? which can be partitioned into cubes. They аге: 


227-10648 with the partition 1, 0, 64, 8, 
and 
303°=27818127 with the partition 27, 8, 1, 8, 1, 27. 


In spite of this the sequence has infinitely many members with the last digit 
#0. This is due to the fact that all cubes of the form (3-107* 2+3) can be 
written in the form 27-10%**+81-10%+81-10%*?+27 which gives the 
partition: 

27, 0, 0 ...(3k zeros), 81, 0 ,0 ...(3k zeros), 81, 0 ,0 ... (3kzeros), 27. 


The integer 303 can be looked at as the generator of a whole family of 
members of the cube-partial digital sequence. 


Question: Is 303 the only generator that guaranties that there are infinitely 
many members of this sequence with non-zero last digits? 


15 Partition of {1, 2, 3, ... 2m} into trigrades using Pythagorean triples. 


The relation 2°+3*+7*=1"+5*+6* holds for k=1,2. It is called а bigrade and is 


written 2317 1,5,6. To find similar relations which hold for a larger range of 
values for k has attracted a lot of interest [4]. Here is an example of a 
pentagrade 


1510182327} 2313152526 


Definition: A relation 

аза ... a, = bi +b, + ... b," 
which holds for m=1, 2, .... k is called a k-grade. Іп abbreviated form this is 
written 

ал, 3,,...а,5 Dy, δν... ba 


The purpose of this study is to give an explicit expression for the partition of 
S={1,2, .. 2n) into a-( a), 2... а.) and B={ b, b; ,...,b,) so that AZ B '. We 
will see that this partition can be obtained through a suitable choice of 
Pythagorean numbers. The trigrades generated in this way will in general 
have a large number of terms. They can, however, be used to produce 
trigrades with only eight terms: 

21,22,23 a, = 5,0% 


Let A, consist of all odd integers x in the interval 1<х<п-1, where n is 
assumed to be even, and all even integers in the interval n*2xx«2n. Во-Ао 
Consider 
2n 
S. -5 ук" (Half the sum of m™ powers of S) 
k=1 
n/2 n/2 
Aa = LQk-1) + £E(n42k)" 
k=} k=1 
Evaluating these sums result in 


! S= AuB with A n B = O, b = а’ (a’ stands for the complement of a). 
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1 
Αι - $ = 0 (i.e. Ao = Bo) (1) 
A2 - S; = n2⁄2 (2) 
Аз - S; = 3n2(2n+1)⁄4 | (3) 
А solution is attempted by constructing new partitions А and B through 
exchange of elements between Ао and Во so that (1) is maintained while 
satisfying (2) and (3). 


Exchange: (xi , x2} = Во — A and {yı , y2} = Ао > B so that 


yi + y> = х + x° (1.1) 
(уу رن(‎ = 7/2 (2.1) 
(y + yz) - (X? + x) = 3п(2п+1)/4 (3.1) 


Assume γ22Υι and write (1.1) as: 

xı نو کس‎ +d; х= у, - d (1.2) 
Substituting this and simplifying (2.1) and (3.1) results іп: 

2d(y; - γι -d) = п2/2 ( 2.2) 

y +у;= 2п +1 (3.2) 
Assume j £ (1, 2, ... 1/2) and write (3.2) as 

y;-n*2j y7n*l-2j (3.3) 
Insert these expressions in (2.2) 


n? - 4d(4j - d -1) (2.3) 
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Integer solutions for d require the discriminant of this equation to be a square 
22. 


(4j - 1 - n” = Z° with n/4 + 1 <j < n/2 (4) 


Let a, b and c be a Pythagorean triple with 2 |b and let k be an odd integer, 
(4) is then satisfied by 4j-l=kc, n-kb and z=ka. It remains to evaluate χι, X2, 
уі and y; in terms of a, b, с and k. To begin with we have j=(kc+1)/4. Since j 
е Z, and ای‎ (mod 4) we can write k=4e-1 where e є Z.. This gives d=(4e- 
1c + a) and the final solution: 


х12 = ((4e - 1X2b + а) + 2 (5) 
yi2 = (de - 120 + c) + 12 (6) 


The conditions 2b>a and 2b>c restrict the number of Pythagorean triples 
which can be used to generate trigrade partitions of the first N positive 
integers where N ts given by: 


N = 2b(4e - 1), еғ Z, (7) 


A further condition is that x, and x2 given by (5) must belong to Ao, i.e. x=0 
(mod 2) if 0<x<n and x=1 (mod 2) if n<x<2n. This is met if 2b+a=-1 (mod 4). 
The corresponding condition must be met for y, and v;. 

From the above conditions and (7) we see that a Pythagorean triple generates 
zero or infintely many trigrades. If it produces infinitely many trigrades then 
the number of terms N of these are in arithmetic progression with first term 
6b and difference 8b. The first non-trivial Pythagorean triples are listed in 
table 9. Those which can be used to produce trigrades are marked in column 
N by giving the number of terms in the trigrade which has the least number of 
terms. 


The trigrade with the smallest number of terms generated in this way is (the 
exchange terms are underscored); 


1, 2, 7, 8, 9, 11, 14, 16, 17, 18, 22, 24 22, 4, 5. 6, 10, 12, 13, 15, 19, 20, 21, 23 


Table 9. Pythagorean triples and their corresponding smallest М. 





Since the sets Αρ and Во are defined from the outset the trigrade partitions can 
be unambigously described in terms of these sets and the exchange terms as 
Shown in table 10. From this table it is seen that the trigrades corresponding 
to the last two Pythagorean triples have the same number of terms. This 
means that we can obtain trigrades with only eight terms by subtraction. For 
example we see from table 10 that 


197, 164, 344, 17 Ы 317, 44, 272, 89 (8) 


А k-grade a;,2;,...a, : bi,b;,...b, is invariant under translation. To prove this 
consider the sets A={xta,, x+a),...,.xta,} and B={x+b, x+b,,...,x+b,}. 
Expand each term to the power m where m<k. 


Хоча” ж Σ Σ jr - B jh jo 2 озь. => 5 9 ж Y xa ы) 
j i=l 


ігі іші i=1j= 


where the assumed k-grade property has been used in exchanging а; and Ы. 
We can consequently reduce all terms in the above trigrade by 17 (if we 
reduced by 16 instead the eight terms would be positive and realively prime 
but would prevent further reductions) to obtain 


0, 147, 180, 327 2 27, 72, 255, 300 (8.1) 
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Тһе k-grade property is evidently invariant under division by а common 
factor. Hence 


0, 49, 60, 109 7 9 24, 85, 100 (8.2) 


which after reduction proved to be identical. 


Table 10. Description of trigrades 


+ ке а =a В 

1 3,4,5 a{1,3,—,11}U{14,16,—,24}U(8,17}-{5, 20} 
24 B={2,4,~,12}U{13,15,—, 23}U{5, 20}-{8,17} 

2 Α-{1,3,..,27}ω{30,32,..,56}ω{18,39}-[{11, 46} 
56 B={2,4,~,28}U{29, 31,..,55}U{11,46}-{18, 39} 

3 ve A={1,3,--,43}( 46, 48,_, 88}U{ 28, 61} - (17,72) 


B-(2,4,—,44)0(45,47,..,87))0(17,72]) - (28,61) 
4 7,24,25 A-[(1,3,.,71)0(74,76,.,144)0(62,83] -(35,110] 


144 B-(2,4, ,72])0(73,75,..,143) (35,110) -(62, 83) 
5 À-11,3,—,167]0(170,172,.., 336) 0(144,193]-(81,256) 
336 B-(2,4,—,168)0(169,171,..,335)0(81,256] -(144,193) 
6 À-(1,3,—,263)0(266,268,..,528)0(226, 303) - (127, 402) 


528. B-(2,4,..,264)0(265,267,.., 527) 0(127, 402] -(226, 303} 


288 — B-(2,4,.,144)0(145,147,..,287)0(35,254])-(62, 227} 

8 À-(1,3,.,335)0(338,340,.,672))2(144, 529) - (81, 592} 
672 B-(2,4,.,336])01337,339,.,671])0(81,592]-(144,529] 

9 A={1,3,-~-,527 }W{530, 532,..,1056])0(226,831) -(127, 930) 


1056 B-(2,4,.,528)0(529,531,.,1055))0(127,930]) - (226,831) 


Table 10. Description of trnigrades, continued. 


+ dm а ma B 


10 11,60,61 A-(1,3,.,179]0(182,184,.,360)0(164,197)-(89,272] 
360 B-(2,4,.,180)0(181,183,.., 359) 0189,272) -(164,197) 
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11 A={1,3,..,419}U{422,424,..,840}U{ 382, 459}-{207, 634} 
| 840 B=(2,4,..,420}U(421, 423,.., 839}U{207, 634}-{382, 459} 
12 91,60,109 A-(1,3,.,179]0(182,184,.., 360}U{44, 317}-{17, 344) 
360 B={2,4,..,180}U{181,183,..,3590(17,344}-(44, 317} 
12 A={1,3,..,419}U{422,424,...,840}U{102,739}-{39, 802} 
‚ 840 B-(2,4,..,420)0(421,423,..,839)0(39,802] -(102, 739} 


Table 11. 65 


ai a2 a3 a4 bi b; b3 Б, 
0 49 60 109 9 24 85 100 
0 67 102 169 22 25 144 147 
0 85 136 221 25 36 185 196 
0 103 122 228 22 45 180 203 
0 107 242 349 25 62 287 324 
0 125 638 763 50 63 700 713 
0 129 244 373 48 49 324 325 
0 151 382 533 49 74 459 484 
0 155 230 385 34 77 308 351 
0 165 184 349 25 84 265 324 
0 181 384 565 60 81 484 505 
0 185 262 447 10 147 300 437 
0 203 554 757 81 86 671 676 
0 215 370 585 46 117 468 539 
0 229 304 533 49 108 425 484 
0 233 278 511 63 86 425 448 
0 241 556 797 72 121 676 725 
0 271 678 949 102 117 832 847 
0 283 542 825 58 165 660 767 
0 293 370 663 75 118 545 588 
0 301 456 757 81 132 625 676 
0 305 458 763 63 158 605 700 
0 331 994 1325 110 169 1156 1215 
0 343 1074 1417 117 174 1243 1300 


аі 82 аз а; bi b> b3 b, 
0 349 376 725 49 180 245 676 
0 359 506 865 81 170 695 784 
0 365 400 765 76 153 612 689 
0 379 522 901 54 225 676 847 
0 381 640 1021 121 156 865 900 
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385 996 1381 96 225 1156 1285 
407 1262 1669 122 225 1444 1547 
461 1392 1853 153 236 1617 1700 
469 1264 1733 108 289 1444 1625 
477 764 1241 153 188 1053 1088 
501 688 1189 76 289 900 1113 
511 1348 1859 148 275 1584 1711 
541 736 1277 121 252 1025 1156 
551 796 1347 147 236 1111 1200 
593 824 1417 121 296 1121 1296 
631 1138 1769 169 306 1463 1600 
649 964 1613 169 288 1325 1444 
787 838 1625 162 325 1300 1463 
817 1048 1865 245 292 1573 1620 
835 1880 2715 296 363 2352 2419 
845 2522 3367 170 567 2800 3197 
849 1084 1933 169 408 1525 1764 
925 2190 3115 34 847 2268 3081 


ΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟ 


L6 Smarandache Non-Null Squares. 


Question: In how many ways can n be written as a sum of non-null squares, 
disregarding the order of the terms? 
Example: 9=12+12+12+12+12+12+12+12+12 

=12+12+12+12+12+22 

=12+22+22 

=32 
Let us denote the number of ways by n(n). The example shows that 
1n(9)-4.Obviously the non-null squares representation which has the largest 
mumber of terms is a representation by n squares of 1 and the representation 
with the smallest number of terms is a representation where we use the largest 
possible squares. If n is a perfect square then (Jn)? is one-term 
representation. 
Let's assume that we know n(x) for all integers x<n-1. Furthermore we make 
the definition: 

lif nis a square 
BUY ыште 


With these precisions we can write the algorithm for calculation of n(n) in 
the following way: 


n(n) = 1 --δ(π) + ae Z (а -i?), with 11) = 1 


This algorithm was implemented for n<100. As is shown in table 12 this 
sequence, after a slow start, grows very rapidly. 


Table 12. The Non-Nu& Square Representations. 





1 1 1 2 2 2 2 3 4 4 
4 5 7 7 7 9 12 13 13 16 
20 23 23 27 35 41 42 47 61 71 
75 82 104 124 134 146 178 217 237 258 
307 377 419 456 535 651 139 804 333 1126 


1300 1422 1629 1955 2275 2513 2846 3397 3972 4435 
4990 5904 6933 7807 8766 10268 12097 13718 15409 17895 
21087 24076 27076 31248 36736 42214 47568 54636 64017 73924 
835546 95596 111637 129306 146714 167379 194807 226021 257447 293255 
340106 394953 451408 514025 594103 690035 790967 901118 1038450 1205451 


L7 Smarandache Non-Null Cubes. 


In how many ways т(п) сап п be written as a sum of non-null cubes? In this 
case there are much fewer representations. For non-null squares we have 
1(9)-4 whereas for non-null cubes we have 1(9)-2. Calculation of the first 
100 terms follows the same method as for the non-null squares with 
corresponding changes in our algorithms. 


lif nisa cube 
0 otherwise 


9-1 


[Vn 2. 


n) =1+8(n)+ > n(n- Ë), with (1) =1 


Table 13. The Non-Null Cube Representations. 
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Chapter II 


Recursive Integer Sequences 


IL 1 The Non-Arithmetic Progression 


We consider an ascending sequence of positive integers a), a2, ... a, such that 
each element is as small as possible and no t-term arithmetic progression is in 
the sequence. In order to attack the problem of building such sequences we 
need a more operational definition. 


Definition: The t-term non-arithmetic progression is defined as the бес: 
{а, :а; is the smallest integer such that a>a;., and such that for Есі there are at 
most t-1 equal differences a, — a, =a, -q, —.-a, -а, } 


1 5-1 


From this definition we can easily formulate the starting set of a t-term non- 
arithmetic progression: 


{1, 2, 3 ..... t-1, ΠῚ} or (a; : azi for i<t-1 and a=t+1 where t23} 
It may seem clumsy to bother to express these simple definitions in stringent 
terms but it is in fact necessary in order to formulate a computer algorithm to 


generate the terms of these sequences. 


Question: How does the density of a t-term non-arithmetic progression Vary 
with t, i.e. how does the fraction a,/k behave for t23? 


Strategy for building a t-term non-arithmetic progression: Given the 
terms a), a>, ... a we will examine in turn the following candidates for the 
term a: 
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аһ = & + d, d=}, 2, 3, ... 
Our solution is the smallest d for which none of the sets 

{8ι, a2, ... ἃς, a * d, atd-e, atd-2e, ... a +d-(t-1)-e :e2d] 
contains a t-term arithmetic progression. 


We are certain that aka exists because in the worst case we may have to 
continue constructing sets until the term a,*d-(t-1)e is less than 1 in which 
case all possibilities have been tried with no t terms in arithmetic progression. 
The method is illustrated with an example in diagram 1. 


In the computer application of the above method the known terms of a no t- 
term arithmetic progression were stored in an array. The trial terms were in 
each case added to this array. In the example we have for d=1, e=1 the array: 
1,2,3,5,6,8,9,10,11,10,9,8. The terms are arranged in ascending order. 
1,2,3,5,6,8,8,9,9,10,10,11. Three terms 8,9 and 10 are duplicated and 11 
therefore has to be rejected. For 4-3, تی‎ we have 1,2,3,5,6,8,9,10,13,10,7,4 
or in ascending order: 1,2,3,4,5,6,7,8,9,10,10,13 this is acceptable but we 
have to check for all values of e that produce terms which may form a 4-term 
arithmetic progression and as we can see from diagram 1 this happens for 
d-3, e-4, so 13 has to be rejected. However, for 4-5, کی‎ no 4-term 
arithmetic progression is formed and e=6 does not produce terms that need to 
be checked, hence a, = 15. 


1 2 3 4 5 é 7 8 9 (ОП 12 13 14 15 





Diagram 1. To find the 9” term of the 4-term non-arithmetic progression. 
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Routines for ordering ап array іп ascending order and checking for 
duplication of terms were included іп а QBASIC program to implement the 
above strategy. 


apoE 
AFAR 





Diagram 2. ax/k for non-arithmetic progressions with 123, 4, 5, ... 15. Bars are shown 
for k = multiples of 10. 


Results and observations: Calculations were carried out for 3<t<15 to find 
the first 100 terms of each sequence. The first 65 terms and the 100^ term are 
shown in table 1. In diagram 2 the fractions a,/k has been chosen as a 
measure of the density of these sequences. The looser the terms are packed the 
larger is a,/k. In fact for {2100 the value of a,/k = 1 for the first 100 terms. 
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In table 1 there is an interesting leap for t-3 between the 64% and the 65” 
terms in that aç, = 365 and ass = 730. Looking a little closer at such leaps we 
find that: 


Table 1. The 65 first terms of the non-arithmetic progressions for f=3 to 15. 


8 ІЗ 4 [55 6 d7 [8 ғ? Бю Би 12 ۱ا‎ ікі4 5 


οσο ч AO +. WH 
O 00 ч ane WH 
οσο σον AWN 
9 оо ч O. Q لد‎ ων” 


OS OO чоо x Q N —‏ سن وې ځن ههو یا دا 


N Ἢ әд πὸ ου o ټټم‎ туф «οὶ — 
هوپ ېو مان ټارن وس بویا کا دا‎ ль ом ~ 


№ 
^ 
м 
ж 


1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
1] 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 


ы 


“SRBSRRESBEKREBRSRIBL 





GREEBRSBIBUEH 
RSEBRLLESKEAVSSRELBRSVBY 
SXTSSRBRREBKREBYRKARBRISSIVSGOK=S о 
5:5%8 ебеввевзчер 

وه یرن ځنم 5 ها کا ناه غانا یم یاه د Κο‏ با غه څ оь‏ 
D‏ کډ تا 8 ټاو د 8 وو £ ها δια‏ 
SIBREBBRESRBBUSR‏ 
SRESKREBKRSESBVSRREBR‏ 

3839559 59698593 α σσ χΩΝ 
ЗесебсеввезызоуваВ 


م 
кэ‏ 


Table continued on next page. 
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Table l continued. 
ї=10 n t12 1=13 t=14 


ΣΙ ΠΩ ΣΤΙΧΟ ΕΞ 


45 
46 
47 
48 
49 
50 
53 
55 
56 
57 
58 
59 
60 
64 
65 
66 
67 
68 
69 
70 
71 
78 
79 
80 


956560666 


O 





14 =3.5-] 28 =2-14 
41 -3.14-1 82 =2.41 
122 =3-41-1 244 =2:122 
365 =3-122-1 730 =2-365 


Does this chain of regularity continue indefinitely? 
ΕΦ. مس کو ادن‎ Ἴτ. 
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Sometimes it is easier to look at what is missing than to look at what we have. 
Here are some observations on the only excluded integers when forming the 
first 100 terms for 1, 12 13 and 14. 


For t=11: 11, 22, 33, 44, 55, 66, 77, 88, 99 The nth missing integer 
is 11.ῃ 

For t=12: 12, 23, 34, 45, 56, 67, 78, 89, 100 The nth missing integer 
15 11-п+1 

For 13: 13, 26, 39, 52, 65, 78, 91, 104 The nth missing integer 
is 29 

For t= 14: 14, 27, 40, 53, 66, 79, 92, 105 The nth missing integer 
is 13-n+1 


Do these regularities of missing integers continue indefinitely? What about 
similar observations for other values of t? 


IL2 The Prime-Product Sequence 


Definition: The terms of the prime-product sequence are defined through (t, : 
t, = p,J* 1, р, is the nth prime number}, where p,# denotes the product of all 
prime numbers which are less than or equal to р,. 


The sequence begins (3, 7, 31, 211, 2311, 30031, ... Y. In the initial 
definition of this sequence t; was defined to be equal to 2. However, there 
seems to be no reason for this exception. 


Question: How many members of this sequence are prime numbers? 


The question is in the same category as questions like ‘How many prime twins 
are there?, How many Carmichael numbers are there?, etc.' So we may have 
to contend ourselves by finding how frequently we find prime numbers when 
examining a fairly large number of terms of this sequence. 


From the definition it is clear that the smallest prime number which divides ما‎ 
is larger than p,. The terms of this sequence grow rapidly. The prime number 
functions prmdiv(n) and nxtprm(n) built into the Ubasic programming 
language were used to construct a prime factorization program for п<10'°. 
This program was used to factorize the 18 first terms of the sequence. An 
elliptic curve factorization program, ECM.UB, conceived by Y. Kida was 
adapted to generate and factorize further terms up to and including the 49th 
term. The result is shown in table 2. All terms analysed were found to be 
square free. À scatter diagram, Diagram 3, illustrates how many prime factors 
there are in each term . 


The 50th term presented a problem. 1;=126173-п, where n has at least two 
factors. At this point prime factorization begins to be too time consuming and 
after a few more terms the numbers will be too large to handle with the above 
mentioned program. To obtain more information the method of factorizing 
was given up in favor of using Fermat's theorem to eliminate terms which are 
definitely not prime numbers. We recall Fermat's little theorem: 


If p is a prime number and (a, p)-1 then а?! = 1 (mod p). 
а"! = 1 (mod n) is therefore a necessary but not sufficient condition for n to be 
a prime number. If n fills the congruence without being a prime number then 


n is called a pseudo prime to the base a, psp(a). We will proceed to find all 
terms in the sequence which fill the congruence 


а" = Ι(πιοά {) 
for 50 < n s 200. t>oo is а 513 digit number so we need to reduce the powers of 


a to the modulus t, gradually as we go along. For this purpose we write t,-1 to 
the base 2: 


1-1 = Y 6()-2! , Where (К) = (0,1) 


From this we have 
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Tobie 2. Prime factorization of prime-product terms 


Prime number 
Pnrne number 
Prime number 
Prime number 


510511 519 - 97 - 277 
9699691 x347 - 27953 
223092871 317 - 703763 
6469693231 2331 . 571 - 34231 
200560490131 
7420738134811 х181- 60611 - 676421 
304250263527211 =61 . 450451 - 11072701 

13062761331670031- 167.78339888213593 
614889782588491411 = 953. 46727 - 13808181181 
32589158477190044731 = 73 - 139 - 173 - 1856476186030) 

19227 60350154212639071 =277 - 3467 - 105229 - 19026377261 
117288381359406970983271 2223 . 525956867082542470777 
7858321551080267055879091 254730729297 - 143581524529603 
557940830126698960967415391 = 1063 - 303049 - 598841 - 2892214489673 
40729680599249024150621323471 = 2521 . 16156160491570418147806951 

121 76447 67340672907 899084554131 = 22093 - 1503181961 - 9688841 4202798247 
2670645! 568927585 135562401 7992791 = 265739 - 10049880359 64897 329167431269 
23768741 896345550770650537 601358311 = 131 . 1039 - 2719- 
6422589 | 88429437337 18061 4] 


23055679 639 455 184247531021 47331756071 = 2336993 - 13848803 - 

7123743602409 1007473549 

23286236435849 7 3607000633 | 688050736307 1 = 960703. 

242387 46455 3038099079 59 4127301057 

239848238289 25228 172706521 638649225839621 | = 2297 - 9700398839 . 

179365737007 - 6001315443334531 

25663761 1759499941 447959781 534007 1 64839447} = 149 -13203797 - 

30501 264491063137 -4276784365108371 1 

27973499 68 | 7854936 1 782761 6187206780967 4997231 = 334507 - 1290433 - 

648046444234299714623177554034701 

3161005464041760778814520629 154366249327 4686991 = 5122427 - 

2025436786007-304670759 5069 540247157055819 

4014476939333036189094441199026045136645885247731 = 

1543 - 49999 - 552001 . 57900988201093 -162808052999907 3967231 
25896479052627740771371797072411912900610967452631 = 

1951.22993.11723231859473014144932345466415143728266617 





Tabie 2. Prime factorization of prime-product terms, continued. 





N=p#+! and its factors 






881 - 1657 - 32633677 - 160823938621 - 5330099340103 . 176429 1759303233 
?| 56 11001 4646650599 1 90057 509233131 6499 4005736633465320043309! = 
6782799 S9OOSS 28882498681 487 - 1476476861 45442451 39224580493 
14921823509 39279 32005887 57366 1 584} 068547 583863326864530411 = 
87549524399 - 650181 61573521013453 - 2621400768441 34219184937113 
60 [22531953499 1831 17732889023622899 200 1 350685 | 6336235654409 191} = 
23269086799 | 80847 . 968321 3481319911991 636641 54180202427 1084713 
353751 669937 174948406357 67087951 74421 2057570647889977 42242987} x 

1381 - 1867 - 8311930927 . 38893867968570583 - 42440201875440880489 1 13304753 
5766152219975951 6590236300353361 343065653840 1 560606631 985606881} = 

1361 - 2141 14727210560829 - 32267019267402210517 - 613228865630544238382107 
9629474207 359839270569 4621 5901 134429 196419 1 3060621 3075415963491271 x 
205590139 - 53252429177 - 7064576339566763 - 

16658990378732521 938085 1 69535089 6256250780509 59 487 4862046961 68398971] = 
62614127 . 26605801 5609361 1580352333 193927566 | 5852809877 2260689062181793 
2981959277793) 21 4269 | 7245346781 042986892551 12174826003064061 414341 5809 
1 2601 - 1651781 - 8564177 - 358995947 - 1525310189119 -5405328664096618954809 
029861252251 
5397346292805549782720214077673687806275517530364350655459511599582614 
91 2107453 . 5634838141 -891415728096410112334489123965712571636329746284 
































329391 =32999 . 175603474759 - 771 48541513247 -2305961 466437323959598530415 
84242331 6227152033 

19896237639 1 690981 640415251 545285 15360273440272182105821 220397 60954139 
10572271 =21639496447 - 7979] 259059 673329495018 877-1 15230777 1 627975804402 
016210177745361 5909 

391955881 491 631233831 61804554421 175259738677 3361 987 48467804 1 8329079654 
0382737 191=521831-50257723 - 1601684368321 -39081 170243262541027 -23875913 
958369977 1585726531 60969521 

7797722041 683461553249199 10632981 387 668799 67899035509 4509 303247 48685} | 
536164700611 =467 -10723 - 57622771 - 5876645549 9458] 45520867 -486325954430 
626096097 19222040521 4947865503847 

164578355079521038773558101 1435590727981 1 6732266964924941 4629852197255 
73413075187091} 1 =1051 - 2179. 16333 - 43283699 - 7531 1908487 -29281 27 10684839 
4609 659667 286646929 343033404487 2907 384889 

367009731 82733191 6465034565550 1 367323398003 1 295533178261946 2457039 
7331115766721 2931 =138678894681 59 -264647 1 42357 1 660867679 1 59849289670356 
4888 1 000360533429306 | 946803757 2880509 


8331 1209 1 2480434503756284637988 1038241 13467 104086031 46546179777 4807729 
2641632790457335111 = 3187 . 31223 - 1737142793. 11463039340315601 .97310450 
3470446969309 1 13 - 43206785807567 1 8923287 5099 500379 
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7 
6 
5 
4 
3 
2 
1 
9 





Diagram 3. Тһе number of prime factors in the first 49 terms of the prime-product 


sequence. 


This product expression for а 


t, 


~! is used in the following Ubasic program to 


carry out the reduction of a^^! modulus t, Terms for which ©(Е)=0 аге 
ignored in the expansion were the exponents k are contained in the array 
Е%О. The residue modulus t, is stored in F. In the program below the 
reduction is done to base А-7. 


100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 


dim EX(1000) 
MzN-1:$—0 
T=1%=0 

while {M-T}>=0 

inc JXT-2*T 

wend 

dec J%:M=M-T\2:inc Ie: ESO (196) — کل‎ 
if M»0 then goto 120 
ΕΞ] 

for ال‎ to % 

A=7 

for K%=1 to E%1J%) 
A=(AA2)@N 

next 

F-F*A:F-FGN 

nexi 
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This program revealed that there аге at most three terms ما‎ of the sequence іп 
the interval 50<n<200 which could be prime numbers. These are: 


Term #75, N=379#+1, Nis а 154 digit number, 
N-17196201054584064334833405683175430195845756358957425604387711050583216 


5523856261308397965147955578800999455782202456522693290629520826275682227 
5663694111 


N-20404068993016374194542464172774607695659797117423121913227131032339026 
1691759299022444537574104587288429298422716055578188216854906766619853898 
3995862280246598588137613940413837615309610314083466556364674016027975521 

231750135686300363861239066166840623542231 17837423905 10526587257026500302 
6968347932485 26734305801 6341 65948702506367 1 76701 2332980646 1 66635537 169754 
29048751 575597 15041 7381 043934255689۱ 24486029492908966644747931 


І #172.N=1Q21#+1.N s q 4 it й 
М-20832554441869718052627855920402874457258652856889007473404900784018145 
7187286244301915872863160885721486313893793092847430199408859808718870830 
26597753881 31 777260588503833 1 62528205231 11233067921 93540483321 70364563007 
17761 6888535712671 50232508655463442766366 180331 2009807 1 1247645589424056809 
033468323906 7 45795726223468483433625259000887 41 1959 1973239736 134883450319 
13058775358684690576 1 4606627 6875058596 1002361 12260054944287636531 


The last two primes or pseudo primes are remarkable in that they are 
generated by the prime twins 1019 and 1021. 


Summary of results: The number of primes q among the first 200 terms of 
the prime-product sequence is given by 6<4<9. The six confirmed primes are 
terms numero 1, 2, 3, 4, 5 and 11. The three terms which are either primes or 
pseudo primes are terms numero 75, 171 and 172. The latter two are the 
terms 1019#+1 and 1021#+1. 


IL3 The Square-Product Sequence 


Definition: The terms of the square-product sequence are defined through {t, 
: ta = (n) 
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This sequence has а structure which is similar to the prime-product sequence. 
The analysis is therefore carried out almost identically to the one done for the 


sequence. We merely have to state the results and compare 
them. 


The sequence begins {2, 5, 37, 577, 14401, 518401, ... }.As for the prime- 

sequence the question of how many are prime numbers has been 
raised and we may never know. There are similarities between these two 
sequences. There are quite a few primes among the first terms. After that they 
become more and more rare. Complete factorization of the 37 first terms of 
the square-factorial sequence was obtained and has been used in diagram 4 
which should be compared with the corresponding diagram 3 for the prime- 
product sequence. 


ο = мо 4. οι 0 + 





Diagram 4. The number of prime factors in the first 40 terms of the square-factorial 
sequence. 


Diagram 4 is based on table 3 which shows the prime factorization of the 40 
first terms in the square-factorial sequence. The number of factors of each 
term is denoted f. The factorization is not complete for terms numero 38 and 


39. A +-sign in the column for f in table 3 indicates that the last factor is not a 
prime. The terms of this sequence are in general much more time consuming 
to factorize than those of the prime-product sequence which accounts for the 
more limited results in this section. Using the same method as for the prime- 


Tabie 3. Prime factorization of square-factorial terms. 


229442532802560001 z101- 227 1708245569901 

38775788043632640001 

760005445655 1997440001 = 29. 109. 2404319663572286441 

171061 22527243 9942400000 1 = 13441 69-12721 70577304043929 

4377631 3669739 505254400000 I = 1 49- 293800762884 1 577533852349 
1265135465055471 7018521 6000001 = 9049-1 3980942259424! 43240713449 
4099038906779 7 2831 4000998400000 1 ع‎ 1078483531837 103551 35404972973 
1479753045347 481921 354360422400000 | =7 10341 -20831 587 15810409256053586 1261 
53919012181 38992768541 7441 68960000000 1=41-10657- 8681 6017- 

348046955409 -4483247 4984 1 

261028437 1 992958 10926909 1 785 1 13600000001 = 

61-157-272557624725} 705240961 774861 76631513 


126337763604459172488624042399 498240000000 } = 337. 8017. 51 4049836440277 481 - 
90967 4823323537849 

668326769 467589022464821 184293345489 600000001 = 509.1 544837 4429- 

8499 40026045327 47687401741723441 

38495621921 3331! 2769397370021 52967 1 17209600000001 


171707860473207588349837- 455289320701 41 406371 6469396531 758248773 
70359079 63854588237 4689246780656 | 195760321 617199 1 0000000000 1 = 
61. 1733 15461- 359525849. 100636381 1265686901 10069- 117459224951 8207759537897 
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Table 3 continued. 


34177) 4 1871 624056287592837 1043383719 1346761650927 193975961 3935941 7700 
01-193- 3217 866100731 693-394521 434431 45645231 476894644096901 291 197410624286 
816576197 

351814 3 11067739 46895230] 22545281450559 425330223858 | 2620552707 1 52831 0538240000000000 

3618413 

8713 
3819014» 1273552658381 61 455835337350607 1 1274031 9954040039721 559090737 4121359 182397440 






N=(nij2+1 and its factors 










000001 =31 7-373-9030 I $653886808.4889 7 828545563235536086347 48201 1761 63219890777 
16189815715361 

138379035 1762182388 1 86847599250 1 52279701 20131 562363084700690457 559040000000 
000000001273. 57986941373 -3269017431 698277804505207 62862495 1781 7882480631366 
14675401 1637734149869 

1894406899 1 562427 68942779 43633734584709 1094440 | 088750629 552452343985257 60000 
090000000۱ = 1 2740636429788 I - 4910557 1 1943381 28021910109- 302797201 1452403842 
82924307 6981 4272052327 643069 



















0000000000000012233-757-15509 190807 49 1 4281 21 70094885906800637 254241 672273179 
0323638834191 845062531 6785821 6021 

39| 93) 4+ | 4160735733984357 43255481 1027341847802665009 490041 649137701 1603858731 6426506 
240000000000000001 261 -10045457577 41 -679001 282693270665470B45882 1 250934897329 
61534689729372855 102292431 5727303361780) 

40196) 4 166571 7749437497 } 8920876976437 4695648426401 5 1 8406663862032 1 85461 739706282409 
984000000000000000001=89. 701- 187100101 949-570306 1 92879859 1 519563156731 42222 
3621393496539544379 492647728 17 483490202551 13441 















product sequence the terms t, in the interval 40<n<200 which may possible be 
primes were identified. There are only two of them: 


І №=(651)2+1, Ni і š 
6902374028907832895045078197262220379290257695327135803427938010402710065 
246438264965962372444657815 1 41285899657 1 53438534056379295 1822384455130747 
80057000000000000000000000000000001 


І $76 М-|764)2%!1 Ni it : 
3555090270010747854202513135770772648194325666925541647977005250280050084 
1772266884421 39166589065164:3920912930369944999452531006264950775782697850 


7198658011625298409931764786386381 15061 7600000000000000000000000000000000 
0001 


П.4 The Smarandache Prime-Digital Sub-Sequence 


Definition: The prime-digital sub-sequence is the set : 
{M=ao+a;-10+a2-10°+...a,10* :M is prime, all digits ao, ау, 8»...ἃν are primes} 
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The first terms of this sequence аге (2, 3, 5, 7, 23, 37, 53,73,.. y. Sylvester 
Smith [1] conjectured that this sequence is infinite. In this paper we will 
prove that this sequence is in fact infinite. Let’s first calculate some more 
terms of the sequence and at the same time find how many terms there is in 
the sequence in a given interval, say between 10“ and 10” The program 
below is written in Ubasic. One version of the program has been used to 
produce table 4 showing the first 100 terms of the sequence. The output of the 
actual version has been used to produce the calculated part of table 5 which 
we are going to compare with the theoretically estimated part in the same 
table. 


Ubasic program 


point 2 

dim А%46),8%44) 

for I%=1 to 6read ASSI) next 

data 1,4,5,89.0 ‘Digits not allowed stored іп А) 
for =1 to 4reod BRIM) лехі 

data 2,3,5,7 ‘Digits allowed stored in 8%() 

for KR=1 to 7 ‘Calc. for 7 separate intervals 
M=0:N=0 

for E%=1 to 4 ‘Only 2.3.5,7 allowed as first digit 
100 P=BR(ER)*1OAKS:PO=P:S=(BR(E%)+1)* 10^K96:gosub 150 


5855585845 


150 while Р<5 

160 P-nxtprm(P):P$sstr[P) ‘Select prime and convert fo string 
inc N 'Count number of primes 

180 1%-епіРУ):С%-0 'C% wil be set to 1 if P not member 

190 for =2 to L% 

200 for JR=! їо ὁ ‘This loop examines each digit of P 

210 if val(rici(P$.I%.1))=AR(JR) then С%=1 

220 nextnext 

230 if CX€-0 then inc MX ‘if criteria filled count member (m%) 

240 wend 

250 retum 
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Tobie 4. The first 100 terms in the prime-digital sub sequence. 


2 3 5 7 23 37 53 73 223 27 
223 257 277 337 353 373 523 557 577 72 
733 757 773 2237 2273 2333 2357 2377 2557 2753 
2777 3253 3257 3323 3373 3527 3533 3557 3727 3733 
5227 5233 5237 5273 5323 5333 5527 5557 5573 5/37 
7237 7253 7333 7523 7537 7573 7577 7723 7727 7753 
7757 22273 22277 22573 22727 22777 23227 23327 23333 23357 

23537 23557 23753 23773 25237 25253 25357 25373 25523 2553/ 
25577 25733 27253 27277 27337 27527 27733 27737 27773 32233 
32237 32257 32323 32327 32353 32377 32533 32537 32573 33223 


Table 5. Comparison of results. 


0.01261 





Theorem: 

The Smarandache prime-digital sub sequence is infinite. 

Proof: 

We recall the prime counting function π(χ). The number of primes PSX is 
denoted T(x). For sufficiently large values of x the order of magnitude of 


την, χ 
r(x) is given by кымы. مر مر یوو سرو‎ Su 


n(a,b,k) the approximate number of primes іп the interval (b-10*,a-10%). 
Applying the prime number counting theorem we then have: 


b 
à Ф 
log10+ — log10 + , ع‎ 


Potential candidates for members of the prime-digital sub sequence will have 
first digits 2,3,5 ог 7, i.e. for a given k they will be found in the intervals 
(2.10% 4.10"), (5.10΄ 6.10 and (7.10*8-10*) The approximate number of 
primes m(k) in the interval (1051091) which might be members of the 
sequence is therefore: 


n(k)-n(4,2.k)*n(6,5,k) (8,7 k) (2) 


The theoretical estimates of n іп table 5 are calculated using (2) ignoring the 
fact that results may not be all that good for small values of k. 


We will now find an estimate for the number of candidates m(k) which 
qualify as members of the sequence. The final digit of a prime number >5 can 
only be 1,3,7 ог 9. Assuming that these will occur with equal probability only 
half of the candidates will qualify. The first digit is already fixed by our 
selection of intervals. For the remaining k-1 digits we have ten possibilities, 
namely 0,1,2,3,4,5,6,7,8 and 9 of which only 2,3,5 and 7 are good. The 
probability that all k-1 digits are good is therefore (4/10)*!. The probability q 
that a candidate qualifies as a member of the sequence is 


Е k-1 
4-5 τ 3) 


ιν... pn 
is therefore given by m(k)=q-n(k). The estimated values are given іп table 5. 
А comparison between the computer count and the theoretically estimated 
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values shows a very close fit as can be seen from diagram 5 where Юр m is 


plotted against k. 


The prime-digital sub sequence 





Diagram 5. Όσιο m as a function of k. The upper curve conesponds to the 


computer count. 


For large values of k we can ignore the terms 


k 
comparison with log 10 in (1). For large k we therefore have 
(a — b)10* 


and (2) becomes 


loga 


log) 
k 


in 


(1?) 


4-10: 


ΜΈΝ 15610 (2) 
Combining this with (3) we get 
5.2 
و‎ )4( 


From which we see (apply for instance l'Hospital's rule) that πι(])-οο as 
ko. A fortiori the prime-digital sub sequence is infinite. 
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Chapter ПІ 


Non-Recursive Sequences 


111.1 Smarandache Primitive Numbers 


Definition: For a given prime p and positive integer n the Smarandache 
primitive number S,(n) is the smallest positive integer such that S,(n)! is 
divisible by р". 


This sequence is important for the calculation of the Smarandache function 
S(n) which is defined as the smallest integer such that S(n)! is divisible by n. 
We note that S,(n)- S(p ). 


It follows immediately from the definition that the sequence S,(n) for a given 
prime p consists of multiples of p. Furthermore the definition implies that the 
sequence is non-descending, 1.6. 5-(п)< S, (n*1). 
It is evident from the definition that 

S,(n)"np for nsp (1) 


Upper and lower bounds for S, (n) have been established by Pal Grónás [1]: 
(p-1)n+1s S,(n)spn (2) 


Let's asssume S,(n-1)<S,(n) and р1 5,(п) then it follows that S,(n*k-1)- 
S.(n+k-2)= ... = S,(ntl)- S,(n) This observation is expressed in the 
following way in Some Notions and Questions in Number Theory [2] " 
Curious property: this is the sequence of multiples of p, each number being 
repeated as many times as its exponent (of power p) is.” 


Algorithms for the calculation of S,(n) are important as a stepping stone for 
the calculation of the Smarandache function. The lack of a closed formula 
makes this a very interesting topic. Three methods will be described, 
implemented and compared. 

Method 1. 

Let's denote S,(n)}=m. We want to calculate the smallest value of m for which 
1-2-3- ... -m=0 (mod p^). This is carried out in the following Ubasic program: 


10 input "p,nP,N 

20 if N<=P then M=N*P:goto 80 

30 Y=1:Z=PAN 

40 while Ү>0 

50 inc M 

60 Y=Y°M@Z ‘The factorial is reduced mod z each time 
70 wend ‘the loop is executed. 

80 print ٧٨ 

90 end 


This program is effective for small values of n. As we move on to large 
exponents the increasingly large modulus slows down the program. A study of 
the compexity of this short and elegant program, which can equally well be 
used to calculate S(n), has been carried out by S.and T. Tabirca (3). 


Method 2. 

When пер we have S,(n)-np. For п>р our potential solutions are multiples of 
p. The first multiple we need to examine is p^. For a given multiple m of p we 
determine first the largest power of p which is less than m, i.e. we determine 
k so that p'«m«p*"!. We then continue by counting the number of factors p in 
m!. This is given by S = [πυρ]! [πυρ] + [πυρ] + ... + [πυρ]. If S<n then we 
proceed to the next multiple of p and so on until we arrive at 5-0 in which 
case m- S, (n). 


The following UBASIC program has been used to tabulate S, (n) as well as 
the Smarandache function in the Smarandache Function Journal [4]. 


10 input ۲۷١ 
20 if N<=P then M=N*P:goto 130 ‘S,(n)=np in this case 
30 І-Р-1 'i =p-1 so that m starts at p? 
40 while S«N 'in line 60 
50 incl 
60 ΜΞΙΡ ‘The starting point for m=p? 
70 8-0 
80 K-floor(log(M)/log(P)) ‘Determine the largest k for which 
90 for Ј=1 to K 'pk«m 
100 S=S+floor(M/PAJ) 'Count and add the number times 
110 nex 'Di occurs in m! 
120 wend 
130 print P,N.M 
140 end 
Method 3. 
This method is based on a theorem due to C. Dumitrescu and V.Seleacu [5]. 
They prove the following: 


Theorem. If for a given prime p the integer n is expressed in the number base 





Бі, b=1+p, Ьз=1+р+р?, ..., by=1+p+ ... p^, resulting in n=cibi+ оду}... 
+ qb, where the first non-zero coefficient csp and с;<р for i<jsk, then 
$у(пу=ср + cap + ... کی‎ 


In the program below the number base is generated through the recursive 
formula by-p-b +1. 


958256554890 


100 


dim С(100),8(100) 

input "p.n P,N 

if N<=P then M=P*N:goto 150 

В(1)=1:1=1 

while М>-В1(1) 

inc | 

B(I)-B(I-1)*P* 1 ‘Calculate b», bs, etc 
wend 

Κ-|-1 

for i=K to 1 step -] 


110 С(1)=м\В(1) ‘Calculate the coefficients 


120 N-res ο C2: s 

130 M=M+C(I}*PAi ‘Express 5.(п) using these 
140 next ‘coefficients 

150 print M 

160 end 

Results. 


Which of these three programs is the most effective? Since 5-(п)-пр when 
п<р all three programs will be equally efficient for n<p. The programs were 
therefore tested by calculating S,(n) for n=ptl, p*2, .. p*50 for 
p=2,3,5,7,...229 (the first 50 primes). The largest number manipulated in the 
programs is therefore 229°”, which is а 559-digit number. The programs 
were fitted as routines in the same in- and output program. Short integers 
were used whenever possible. To calculate these 2500 values of S (n) method 
I took 31 m 26 s, method II 22 s and method HI 11 s on a Pentium 100 Mhz 
laptop. The time t was measured in milliseconds. In figure 1 In t is plotted 
against the number of primes for which the programs have been executed. It 
is seen that method III is the most effective. It runs about twice as fast as 
method II. Method I is effective for small primes but slows down considerably 
when the modulus р increases. 











өз вз چ‎ азе за 
Prime number numero 


Figure 1. Comparison of execution times for methods I to III. 
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Table 1. Spin} for n=}, 2. ... 48 and р<47 












52 









6 > 3 5 7 H B 7 9 2 9 3 J 41 4 ته‎ 
2 ИЕ & 10 14 2 2 34 3 4 SB 62 74 8 86 94 
3 Ш 15 3 3 8; 93 111 123 129 141 
4 ¿ 20 44 52 116 124 18 164 172 188 
5 δ "E= 55 65 145 155 185 205 215 235 
6 8 15| 25 66 78 174 186 222 246 258 282 
7 з 1з з 7 9 203 217 29 287 301 329 
8 10 18 35 88 104 232 248 296 328 344 37 
9 12 2 4 56 17 261 279 333 349 387 423 
12 24 45 63 130 290 310 370 40 430 С 
14а 27 5 143 39 341 407 451 473 517 
ЕО Ж 59 348 372 444 492 516 564 
27 55 377 403 481 533 559 611 
E 30 ته‎ 406 434 518 574 602 658 
d 33 65 435 465 555 615 645 704 
I 70 464 496 592 656 688 7 
' $ 75 43 527 629 697 731 
20 39 75 522 558 666 738 774 

22 42 50 551 589 703 79 δι; € 

24 45 85 580 620 70 820 860 940 
24 45 9 609 61 777 861 903 987 
24 48 95 638 682 814 902 946 1034 
26 51 667 7۱3 851 943 989 1081 

28 em 696 744 888 984 1032 1128 
28 54 105 456 552 725 775 925 1025 1075 1175 

30 54 0 475 575 754 806 962 1066 1118 1222 
27 32 57 115 494 598 783 837 999 1107 1161 1269 
25 32 60 120 175 286 338 459 513 621 812 848 1036 1148 1204 1314 
2% 32 63 125 182 27 351 476 532 644 899 1073 1189 1247 134 
مد‎ 32 63 125 189 308 364 493 551 667 930 1110 1230 1290 141 
31 32 66 Es 39 377 510 570 690 870 1147 1271 1333 1457 
32 34 69 m 330 390 527 589 713 899 1184 1312 1376 1504 
33 46 72 135 203 341 403 544 66 736 928 992 1221 1353 1419 155 
34 36 72 10 210 352 416 561 627 759 957 1023 1258 1394 1462 1598 
35 38 75 145 217 363 429 578 646 782 986 1054 1295 1435 1505 1645 
36 0 78 10 224 363 442 528 665 805 1015 1085 1332 1476 1548 165 
37 40 81 150 231 374 455 595 684 828 1044 ЛЕ Ы 1517 1591 1739 
38 40 81 155 238 385 468 612 703 851 1073 ШЕ 1558 1634 4 
39 42 81 10 245 396 481 627 722 874 1102 1178 1406 1599 1677 183 
40 44 81 165 245 407 494 646 722 897 1131 1209 1443 1640 1720 86 
41 44 84 170 252 418 507 663 741 90 1160 1240 1480 1927 
42 46 87 175 259 49 507 60 760 943 1189 1271 1517 1974 
43 48 90 175 266 440 520 697 779 966 1218 1302 1554 2 2021 
44 48 90 180 273 451 533 714 798 989 1247 1333 1591 176 2068 
45 48 93 185 280 462 546 731 817 1012 1276 1364 1628 1804 1892 1 
46 48 96 190 287 473 559 748 836 1035 1305 1395 1665 1845 1935 216 
ai so 99 195 29 572 765 855 1058 1334 1426 1702 1886 19 (ЕН 
48 52 99 20 29 84| 585 782 874 1058 1363 1457 1739 1927 2021] 2209 


In table 1 consecutive squares have been emphasized by frames. The series іп 
light frames follow the rule S,(p)À- S,(p+1)=p* which is obvious from the 
definition. The series in bold frames is less obvious and this study will be 
concluded by proving the following: 
S,(4p+4y= S(4p*3)-4p' or δρ δρ ΩΡ} 
Proof: Applying the theorem of Dumitrescu and Seleacu we have: 
4p+4—4-(1+p), i.e. c,=0 and «7-4 from which 


S,(4ptp)=0-pt4-p=(2p)” 


4p+3=p-1+3-(1+p), i.e. ci=p and «732 from which 


S (4p*3)-p-p*3-p - QpY 


III.2 The Smarandache Function S(n) 

Definition: S(n) is the smallest integer such that S(n)! is divisible by n. 

The properties of this function have been subject to detailed studies in many 
papers. Some of these properties are listed here. When needed n is 
represented in the form n = p! پو دم‎ 


1. (msn 
2. pis prime iff S(p)-p 


3. If m-nx0, then S(mn)sS(n)-S(m) 
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4. If (m.ny=1, then 5(шп)-шах(5(ш),5(п)) 
5. S(m=max(S, (αι), S, (%2), ---Sp, (x )} , where S,(n) are 
Smarandache primitive numbers as defined in section П.1 


6. Hn is square free, Le. a,-a,= ...7o&, then Ә(пу-шахір,, po, -..Ρι) 


The purpose of this section, however, is to give an effective computer 
algorithm! for the calculation of S(n) and to provide a table for the first 1000 
values of this function. 


10 ‘A UBASIC program for online caiculation of S(n) 
20 dim Е(10),Е%(20) 

30 input “N "N 

40 gosub 70 

50 print "S(n)s “S 

60 end 

70 NI=N:gosub 170 


N is factorized and the prime factors stored іп Ғ/) with the 
corresponding exponents in ΕΘ} 


80 if F{1}=N then S=N:goto 160 


If F(1)=1 then n is prime and consequently S(nj=n 

90 for J%=1 to КФ 

100 if E%(J%)}>1 then P2=F(J%):N2R=ER(J%):gosub 260:F[J%)=M 
110 next 


For prime factors of n with exponents larger than 1 F{J%]=P2 is 
replaced by ҒҚ(Ж/ғМ which Б calculated Бу one of the routines 
described есте (here method 2 J 


120 5-0 


1 The author was inspired to do this after several requests {ο heip provide values for this function for 
large values of n. 


130 
140 
150 


for JR=1 to K% 
if F(J%)>S then S=F(J%) 
next 


The largest value of F[] is the valve of S{n} 


160 
170 


retum 
'Factorisation of N1, output f[i) Ae% (i) 


A simple routine fo іосіопге Ni 


180 
190 


K%=0:for I%=1 to 10:F(I95) -0:next 
for 125-1 to 20:E9e(175) 20: next 
while N1»1 
-prmmdiv(N1) 
if P=F(K%) then inc EX(K%) else inc КЖ%:Ғ(К%)-Р:іпс EX(K%) 
سوا وب د‎ 
wend 
return 
‘Calculation of Sp(n), In: p2,n2%, Out: m. This routine 6 
documented in section ή, 1 
12=0:52%=0 
while 52%<М2% 
іпс 12 


K2=floor(log(M} /log(P2}) 
for 12-1 to K2 
S2%=S2%+floor[M/P2AJ2) 
next 

wend 

return 


This is the bare minimum program. When used in a special problem the 
number of instructions tend to increase considerably because of “safety 
valves” and data flow from and to other processes. 
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Table 2 ctd. Тһе Smarandache function 





n : 0 1 2 3 4 5 6 7 8 9 
51: 17 73 12 19 257 103 43 47 37 173 
52 : 413 521 29 523 131 10 263 31 11 46 
53: 53 59 19 41 89 107 67 179 269 14 
54 : و‎ 541 27 181 17 109 13. 547 7 61 
55: 3 29 23 79 277 37 139 557 31 43 
se  : 7 17 281 563 47 113 283 9 71 569 
57 і 19 571 13 191 41 23 8 57 34 193 
58: 29 83 7و‎ 53 73 13 293 587 14 31 
55. : 59 197 37 593 11 17 149 199 22 599 
ФГ: 10 601 43 67 151 22 101 607 19 29 
б: 61 47 17 613 307 41 11 617 103 619 
6: : 31 22 31 89 13 20 313 19 15? 37 
es ; 7 631 29. 241 317 7 53 14 29 71 
64 ; 8 641 107 3 23 43 19 647 9 59 
6  : 13 31 163 653 109 131 41 73 47 659 
6 : 11 661 331 17 83 19 37 29 167 223 
€7 i 6 61 8 673 337 10 26 677 113 7و‎ 
68 : 17 227 31 683 19 137 21 229 43 53 
69 : 23 691 173 11 347 139 29 41 349 233 
70: 10 701 13 37 1 47 353 101 59 709 
72: 71 79 89 31 17 1 179 239 359 719 
72: 6 103 38 241 181 29 22: 727 13 15 
7323: 73 43 61 733 367 14 23 67 41 739 
74, i 7 19 583 3 31 149 373 83 17 107 
725 |І: 15 3 47 251 29 151 9 757 379 23 
746 1 19 761 127 109 191 17 383 59 10 769 
277: 11 257 193 773 43 31 97 37 399 41 
794: 13 71 23 29 14 157 131 787 197 263 
359. 1 79 113 21 61 397 53. 199. 797 19 47 
80: 10 89 401 73 67 23 31 269 101 809 
οι  : 9 811 29 271 37 163 17 43 409 13 
82 : 41 821 137 823 103 11 59 827 22 829 
83 : 83 277 13 17 139 167 19 31 419 839 
94  : 7 58 421 281 211 26 47 22 53 283 
85 і 17 37 71 853 61 19 107 7 13 859 
$6 : 43 41 431 3 9 173 433 34 31 79 
87: 29 67 109 97 23 25 з 877 439 293 
88 : 11 881 14 883 17 59 443 887 37. 127 
89 : 89 11 223 47 149 179 8 23 449 31 
90 ; 10 53 41 43 113 181 151 907 227 101 
э1_ : 13 911 19 83 457 61 229 131 17 919 
92 : 23 307 461 ΤΙ 11 37 463 103 29 929 
эз : 31 19 233 311 467 17 13: 937 67 313 
94 : 47 91 157 41 59 9 43 947 79 73 
95 : 19 7 17 3 53 191 239 29 479 137 
96 {| 8 62 37 107 241 193 22 967 22 19 
97 : 97 97 12 139 487 13 61 977 163 89 
98 : 14 109 491 983 41 197 29 47 19 43 
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37 499 997 83 199 71 331 31 991 11 : وو 
Smarandache m-Power Residues‏ 111.3 


Definition: The m-power residue of n is the largest m-power free number 
which divides n. Notation Μη). 

Let's express n in the form n = рг! -ρ)7-...Ρν' then M, (n) = p? різ. ραν ; 
where о; -min(a; ,m-1) for i=1, 2, ... К. 


It follows directly from the definition that М.(п) is а multiplicative arithmetic 
function, which is expressed in (1) below. 


H (n,m)-1 then M.(nm)-M.(n)-M.(m) 


We will only make calculations for two special cases: Cubical residues and 
Square residues. 


IIL3a Cubical Residues. 


In the prime factor representation of n تجه‎ will be replaced by о;-2. А 
simple Ubasic routine results in table 3. 


Tabie 3. The first 100 terms of the cubical residues sequence. 








n : 0 1 2 3 4 5 6 7 8 9 

; 1 2 3 4 5 6 7 4 5 
1 і 10 11 12 13 14 15 4 17 18 19 
2 1 20 21 22 23 12 25 26 9 28 29 
3 | 30 31 4 33 34 35 36 37 38 39 
4 |І 20 41 42 43 44 45 46 47 12 49 
6 i 50 51 52 53 18 55 28 57 58 59 
6 : 60 61 62 63 4 65 66 67 68 69 
7: 70 71 36 73 74 75 76 77 78 79 
8 } 20 9 82 83 84 85 86 87 44 89 
9 |І 90 91 92 93 94 95 12 97 98 99 


TIL3b Square Residues 


For n in the form n=p;'-p?-...py* the square residue S,(n) is given by 
Sí(Dn)-prpz ... р. 
Apart from the multiplictive property we add two more properties of this 
sequence. 
If a>0 and p any prime number then و( م5‎ 
S(am)=S,(n)-S,(m) /(S.(n), S,(m)). 
The first 100 terms of this sequence are given in table 4. 


Table 4. The first 100 terms of the square residues sequence 








n 0 1 2 3 4 5 6 7 8 9 

1 2 3 2 5 6 7 2 3 
1 10 11 6 13 14 15 2 17 6 19 
2. 10 21 22 23 6 5 26 3 14 29 
3 30 31 2 33 34 35 6 37 38 39 
4 10 41 42 43 22 15 46 47 6 7 
5 10 51 26 53 6 55 14 57 58 59 
6 320 61 62 21 2 65 66 67 34 69 
7. 70 71 6 73 74 15 38 77 78 79 
8 2 3 82 83 42 85 86 87 22 89 
9 30 91 46 93 94 95 6 97 14 33 





In chapter VI we will return to the Smarandache function and the square 
residues. 
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Chapter IV 


Periodic Sequences 


IV.1 Introduction 


In Mathematical Spectrum, vol 29 No 1 [1], is an article on Smarandache's 
periodic sequences which terminates with the statement: 


*There will always be a periodic sequence whenever we have a 
function f:S—S, where S is a finite set of positive integers and we 
repeat the function f. " 


However, consider the following trivial function f(x,):S S, where S is an 
ascending set of integers (ai, a2, ... a, ... an: 


Í xi ifx>a, 
(х=! x, Ғхұға, 
( ха Ғхұ<а, 


As we сап see the iteration of the function f in this case converges to ап 
invariant a,, which we may of course consider as a sequence (or loop) of only 
one member. In this study, however, a periodic sequence will be referred to 
as such if it has at least two members. If it has only one member it will be 
referred to as invariant. 


There is one more snag to overcome. In the Smarandache sequences 05 15 
considered as a two-digit integer. The consequence of this is that 00056 is 
considered as as a five digit integer while 056 is considered as a three-digit 
integer. We will abolish this ambiguity, 05 is a one-digit integer and 00200 is 
a three-digit integer. 


With these two remarks in mind let’s look at these sequences. There аге in all 
four different ones reported in the above mentioned article in Mathematical 
Spectrum. The study of the first one will be carried out in much detail in view 
of the above remarks. 


IV.2 The Two-Digit Smarandache Periodic Sequence 

It has been assumed that the definition given below leads to a repetition 
according to Dirichlet’s box principle (or the statement made above). 
However, as we will see, this definition leads to a collapse of the sequence. 


Preliminary definition. Let N, be an integer of at most two digits and let N,’ 
be its digital reverse. We define the element Мұ, of the sequence through 


ΝκηΞ М-М | 
where the sequence is initiated by an arbitrary two digit integer №. 
Let's write N, in the form N,=10a+b where а and b are digits. We then have 
N= | 10a+b-10b-a | =9. | а | 
The |а | сап only assume 10 different values 0,1,2, ... ,9. This means that 
N; is generated from only 10 different values of Ν.. Let's first find out which 
two digit integers result in [а |= 0,1,2, .. and 9 respectively. 


|αοἱ 554 two digit integers 


1] 77 99 

10 2 i A - 5 43 i: 54 56 65 67 76 78 87 89 98 
13 20 24 31 35 42 46 53 57 64 68 75 79 86 97 

14 25 30 36 41 47 52 58 63 69 74 85 96 

۱5 26 37 40 48 51 59 62 73 84 95 


27 38 49 SO 61 72 
7 2 39 60 71 82 93 
19 2 70 81 92 

19 80 91 


ON о ο Á WH — 0‏ نلا 
ο.‏ 
ο‏ 
^ 
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It is now easy to follow the iteration of the sequence which invariably 
terminates in 0, table 1. 
Table 1. eration of sequence according іо the preliminary definition 


N2 Ns N4 Ns Ns Ns‏ امو 
0 0 
0 9 ] 
0 9 45 27 63 18 2 
0 9 45 27 3 
0 9 45 27 36 4 
0 9 45 5 
0 9 54 6 
0 9 45 27 63 7 
0 9 45 72 8 
0 9 45 27 63 81 9 


The termination of the sequence is preceded by the one digit element 9 whose 
reverse is 9. The following definition is therefore proposed. 


Definition of Smarandache’s two-digit periodic sequence. Let N, be an 
integer of at most two digits. №” is defined through 


| (the reverse of N, if N, is а two digit integer 
š: y C RR 
We define the element М of the sequence through 
Νε” NENE | 


where the sequence is initiated by an arbitrary two digit integer N, with 
unequal digits. 


Modifying table 1 according to the above definition results in table 2. 
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Conclusion: The iteration always produces а loop of length 5 which starts оп 
the second or the third term of the sequence. The period is 9, 81, 63, 27, 45 or 
a cyclic permutation thereof. 


Table 2. Iteration of the Smarandache two digit sequence 


а-ы N2 N3 М, М5 М5 Ne № 
1 9 81 63 27 45 9 
2 18 63 27 45 9 81 63 
3 27 45 9 81 63 27 
4 36 27 45 9 81 63 27 
5 45 9 81 63 27 45 
6 54 9 81 63 27 45 9 
7 63 27 45 9 81 63 
8 72 45 9 81 63 27 45 
9 81 63 27 45 9 81 





IV.3 Тһе Smarandache n-digit periodic sequence. 


Let’s extend the definition of the two-digit periodic sequence in the following 
way. 


Definition of the Smarandache n-digit periodic sequence. 


Let Мұ be an integer of at most n digits and let R, be its reverse. Мұ” is 
defined through 


N, = R, . 1077! Hes Νε] 
We define the element Мұ, of the sequence through 
ΝΜΗΞ || 
where the sequence is initiated by an arbitrary n-digit integer N, in the 


domain 10"€«N,«10"". It is obvious from the definition that 0<N,<10™", 
which is the range of the iterating function. 
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Let’s consider the cases n=3, п-4, n=5 and п-б. 
n-3. 


Domain 100<N,<999. . The iteration will lead to an invariant or a loop 
(periodic sequence) There are 90 symmetric integers in the domain, 101, 
111, 121, ...202, 212, .., for which М;-0 (invariant) All other initial 
integers iterate into various entry points of the same periodic sequence. The 
number of numbers in the domain resulting in each entry of the loop is 
denoted s in table 3. 


Tabie 3. Smarandache 3-digit periodic sequence 


5 239 11 200 240 120 
Loop 99 891 693 297 495 


It is easy to explain the relation between this loop and the loop found for n-2. 
Consider N=ao+10a,+100a,. From this we have ٨-٧ 99 هه‎ 11-91a;-aol 
which is 11 times the corresponding expression for n-2 and as we can see this 
produces a 9 as middle (or first) digit in the sequence for n=3. 


B=. 


Domain 1000<N,<9999. The largest number of iterations carried out in order 
to reach the first member of the loop is 18 and it happened for N,=1019. The 
iteration process ended up іп the invariant 0 for 182 values of М, 90 of these 
are simply the symmetric integers in the domain like N,=4334, 1881, 7777, 
etc., the other 92 are due to symmetric integers obtained after a couple of 
iterations. Iterations of the other 8818 integers in the domain result in one of 
the following 4 loops or a cyclic permutation of one of these. The number of 
numbers in the domain resulting in cach entry of the loops is denoted s in 
table 4. 


] This is elaborated in detail in Surfing on the Ocean of Numbers by the author, Vail Univ. Press 
1997. 


Table 4. Smarandache 4-digit periodic sequences 


5 378 259 
Loop 2178 6534 
5 324 18 288 2430 310 
Loop 90 810 630 270 450 
5 446 2 449 333 208 
Loop 909 8181 6363 2727 4545 
5 329 11 290 2432 311 
Loop 999 8991 6993 2997 4995 


ηΞ 5, 


Domain 10000<М,<99999. There are 900 symmetric integers іп the domain. 
920 integers in the domain iterate into the invariant 0 due to symmetries. 


Table 5. Smarandache 5-digit periodic sequences 


5 3780 2590 
Loop 21978 65934 

5 3240 180 2880 24300 3100 

Loop 990 8910 6930 2970 4950 

5 4469 1] 4490 3330 2080 

Loop 9009 81081 63063 27027 45045 

5 3299 101 2900 24320 3110 

Loop 9999 89991 69993 29997 49995 

n=6. 


Domain 100000<N;<999999. There аге 900 symmetric integers іп the 
domain. 12767 integers in the domain iterate into the invariant 0 due to 
symmetries. The longest sequence of iterations before arriving at the first loop 
member is 53 for N=100720. The last column in table 6 shows the number of 
integers iterating into each loop. 
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Tabie 6. Smarandache 6-digit periodic sequences 


| 
11 0 

L2] 13586 6530 

| 
ι 3 J 219978 659934 

L 900 8100 6300 2700 4500 


4 

5 25057 18 12259 20993 4449 62776 
L5 9090 81810 63630 27270 45450 

5 47931 14799 42603 22941 29995 
L6 9990 89910 69930 29970 49950 

5 25375 11 32375 21266 4405 ( 
L7 | 90009 81008! 630063 270027 450045 


5 1488 2 1005 1033 237 37 
18) 90909 818181 636363 272727 454545 


1809 11 1350 1570 510 250 
19 | 99099 891891 693693 297297 495495 


Ж 1919 2648 7292 123573 12472 
[10] 99999 899991 699993 299997 499995 


5 
1 
152 4 1254 972 492 111 826 485 429| 4725 
Ln 
5 
1 











M Ze 


Ë 


10989 978021 857142 615384 131868 736263 373626 252747 494505 
623 64 156 796 377 36 525 140 194 
43659 912681 726462 461835 76329 847341 703593 308286 374517 
596 117 156 793 327 65 530 139 179 
340956 318087 462726 164538 670923 341847 406296 286308 517374 



















L 12 
ctd 
ctd 


Table 7. The subtraction periodic sequence, 10<N:<99 





31 23 42 34 53 45 64 56 75 67 86 78 97 89 9 1 1220 #1 
41 24 52 35 63 46 74 57 85 68 96 79 8 90 19 2 30 13 #2 
51 25 62 35 73 47 84 58 95 69 7 80 18 91 29 3 40 14 #3 
61 26 72 37 83 48 94 59 6 70 17 81 28 92 39 4 50 15 %4 
71 27 82 38 93 49 5 60 16 که 





IV.4 Тһе Smarandache Subtraction Periodic Sequence 


Definition: Let М, be a positive integer of at most n digits and let Б, be its 
digital reverse. М,” is defined through 
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М, = R, 107 Dee ^] 
We define the element Nin of the sequence through 

N= |N اع‎ 
where c is a positive integer. The sequence is initiated by an arbitrary positive 
n-digit integer №. It is obvious from the definition that 0<N,<10""', which is 
the range of the iterating function. 
c1, т-2,10<М,<99 
When N; is of the form 11-k غه‎ k-1 then the iteration process results in 0, 
see figure la. 


Every other member of the interval 10<N, وود‎ is а entry point into one of five 
different cyclic periodic sequences. Four of these are of length 18 and one of 
length 9 as shown in table 7 and illustrated in figures 10 and ic, where 
important features of the iteration chains are shown. 


33 37 99 
82 72 98 
27 26 سه 88 

«- (= 87 4 71 
77 15 16 

> 76 50 60 

66 -ί-]) > 04 «- > 05 «- 
65 —(-) 39 (+9) (e 49 (+9) 
55 92 93 
54 28 38 
44 81 
43 17 
33 70 
32 > 0% -- 
22 (e 59 —>{+9} 
21 94 
1 48 
10 83 
0 37 


Fig. la Fig. Ib Fig. Ic 
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1<с<9, n-2, 100<М,<999 


A computer analysis revealed a number of interesting facts concerning the 
application of the iterative function. 


There are no periodic sequences for c=1, c=2 and c=5. All iterations result in 
the invariant 0 after, sometimes, a large number of iterations. 


Table 8. Loop statistics, L=length of loop, f=first term of loop 

















5 ft/la 0 " 22 33 50 100 167 189 200 
1 N: 900 
2 N: 900 
3 М; 241 59 150 
1 240 
2 210 
4 N: 494 42 
1 364 
5 N: 900 
6 м 300 59 84 
1 288 
2 169 
7 М: 109 535 
1 101 
2 10] 
3 14 
4 14 
5 13 
6 13 
8 Nr 203 43 85 
1 252 
2 305 
3 12 
9 N: 21 79 237 170 
4 20 
5 10 
6 161 
7 12] 
8 81 


For the other values of с there are always some values of N, which do not 
produce periodic sequences but terminate on 0 instead. Those values of Ni 
which produce periodic sequences will either have N, as the first term of the 
sequence or one of the values f determined by 1<f<c-1 as first term. There are 
only eight different possible value for the length of the loops, namely 11, 22, 
33, 50, 100, 167, 189, 200. Table 8 shows how many of the 900 initiating 
integers in the interval 100<N,<999 result in each type of loop or invariant 0 
for each value of c. 


۸ few examples: 
For c=2 and N,=202 the sequence ends іп the invariant 0 after only 2 
iterations: 

202 200 0 
For с=9 and Νι-208 a loop is closed after only 11 iterations: 

208 793 388 874 469 955 550 46 631 127 712 208 
For c=7 and N;7109 we have an example of the longest loop obtained. It has 
200 elements and the loop is closed after 286 iterations: 


109 894 491 187 774 470 67 753 350 46 633 329 916 612 209 895 591 188 874 471 
167 754 450 47 733 330 26 613 309 896 691 189 974 472 267 755 550 48 833 331 
126 614 409 897 791 190 84 473 367 756 650 49 933 332 226 615 509 898 891 191 
184 474 467 757 750 50 43 333 326 616 609 899 991 192 284 475 567 758 850 51 
143 334 426 617 709 900 19 476 667 759 950 52 243 335 526 61 1 


102 194 477 7677 336 626 619 902 195 584 478 867761 
160 54 443 337 726 620 19 903 302 196 684 479 967 762260 55 543 338 826 621 
119 904 402 197 784 480 77 763 360 56 643 339 926 622 219 905 502 198 884 481 
177 7644 7 19 199 984 482 277 76 58 | 
136 624 4199077 493 387 776 670 69953 352 246 635 529 918 812211 
105 494 487 777 770 70 63 353 346 636 629 919 912 212 205 49 587 778870 71 
163 354 446 637 729 9 13 305 496 687 779 970 72 263 355 546 638 829 921 


] 14 405 497 787 780 80 73 363 356 646 639 929 922 222 215 505 498 887 781 


180 74 463 357 746 640 39 923 322 216 605 499 987 782 280 75 563 358 846 641 
139 9244 


24 422 217705 500 2 


IV.5 The Smarandache Multiplication Periodic Sequence 
Definition: Let с>1 be a fixed integer and Νο and arbitrary positive integer. 


Мұн is derived from N, by multiplying each digit x of N, by c retaining only 
the last digit of the product cx to become the corresponding digit of Νιμ. 
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In this case each digit position goes through а separate development without 
interference with the surrounding digits. Let’s as an example consider the 
third digit of а 6-digit integer for تک‎ . The iteration of the third digit follows 
the schema: 

xx7yyy ---- the third digit has been arbitrarily chosen to be 7. 

xxlyyy 

xx3yyy 


XX9yyy 

xx7yyy —-- which closes the loop for the third digit. 
Let's now consider all the digits of a six-digit integer 237456: 

237456 

691258 

873654 

419852 

237456 —— which closes the loop. 
The digits 5 and 0 are invariant under this iteration. All other digits have a 
period of 4 for c=3. 


Conclusion: Integers whose digits are all equal to 5 are invariant under the 
given operation. All other integers iterate into a loop of length 4. 
We have seen that the iteration process for each digit for a given value of c 
completely determines the iteration process for any n-digit integer. It is 
therefore of interest to see these single digit iteration sequences: 


With the help of table 9 it is now easy to characterize the iteration process for 
each value of c. 


Integers composed of the digit 5 result in an invariant after one iteration. 
Apart form this we have for: 


c=2. Four term loops starting on the first or second term. 
c=3. Four term loops starting with the first term. 


c=4. Two term loops starting on the first or second term (could be called a 
switch or pendulum). 


Table 9. One-digit multiplication sequences 


ч CN со Ὃ 
O сч соч O ^O CN OO т 
ч OO N OQ O TWN ^O 


— NM * هب‎ OQ م‎ O ох 


со xO O. CN D ته‎ ~ ч мМ 


O М CO O.‏ زب په CN C)‏ سم 


VON ^O O TON ο 
сч ч ^O 00 O CN ч хо 0 


-- ÇN په زم‎ 10 0 М с O. 


— NM با به‎ OR ته‎ Ον 


oo e CN ч‏ زي 


о сч ооч O ^O (CN O ч 


ORM соо!‏ با * ÇN MM‏ سم 


O L‏ ما O 0 O O‏ زي 
O D»‏ را O VOU O‏ 9 


m NO ч U Ὁ мМ CO O. 


ο.‏ ته NM c `o М‏ سم 


о ΟΡ ^O LO TON — 


O.‏ ته N C) чао ою‏ سم 


TON ^O O TWN Ὃ 


OO о * N O DOO WKN 


O‏ ته ο‏ هی D‏ په ÇN CO)‏ سم 
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Periodic Sequence 


Definition. Let No be a two-digit integer а,-10+ао. If a;+ao<10 then b,= a, +ao 


otherwise b= atat]. bo=la)-aol . We define N, %, -10+bo. New; is derived 


from N, in the same way.” 
There are no invariants in this case. 36, 90, 93 and 99 produce two-element 


loops. The longest loops have 18 elements. A complete list of these periodic 


с=9. Two term loops starting with the first term (pendulum). 
sequences 15 presented below. 


c=7. Four term loop starting with the first term. 
с=8. Four term loop starting with the second term. 


IV.6 The Smarandache Mixed C 


c=5. Invariant after one iteration. 


S 3 3 9 29 
ورس و‎ “κ ч 99 
208924 2 8 
SSR ΝΕΡΟ 
25425 9 92 
2 9 K ç 9 s 2 o 
9030807 2 сч 


Πα 8 8 8 په‎ ١ م‎ 
29098598053 саяс 
ه٣ 7ه"‎ ο5- 8 ὢ 
(80925828 DR OR 


ος προς ορ BABA 


89892392 9R$33929 
88с824258%9сРЕКК98 
383:589<0с2832% 

9:239927229838сғ 
5%942529283%8434 
RQRSSSETSESTESS3I3SS 
- وم8۹‎ ه٥‎ 24595 9 
OQ وه‎ πώ مهه‎ 


pana ыз γρ ұшты y 


2 Formulation conveyed to the author: “Let N be a two-digit number. Add the digits, and add them 
again if the sum is greater than 10. Also take the absolute value of their difference. These are the first 


and second digits of N;.” 


L in 5 K o © 

оз 2 Ὁ I 3 K 83 9 ο 8 E ο 
39 8 59 428 R 94 SF ὦ o 8 0 3943 29 33 8 
8% 2 сю ES JR 05% Sd 9982 a2 S SF oo ОБ 
Ko t 29 ον πα a d 5989 PRR 89% 85 88 oS 
де о no ЧӘ қ 92% Fo d) c 288 Fo Ὦ τ ZS ЖЫН» СЕ 
eas 5 Qo 58 9 K o» олы & 5 8 oo ы оа 2 9 
Quo d mor. ۳٣ S бс SR 3 9 SR SR 3 Qo DD S 
SO к 2R 33 95 FILS Baa 59 299% © 28 95 
om d өз no VN 5 238 RRS 3: 889 BO uot. о 
се و‎ 28 89 v3 %99 CER Бесс БоБо 2 Le Ss 
38 S ЯК RS Яс CSS DAL 8839 RSIS Ὁ 88 So 
ne πό R 25 е oec Baas τς 8948 оз 88 38 
лр с © ӨСЗ onm 4309066995539 Qs د‎ 9 NI JU p2 
ΕΕΣ Ὁ 93 82 DD ο... ος SA 82 
53532 KY 55 98 موک لوان ده لم مه هه ٧ه د ته 6م‎ 08٨ AN 33 وم‎ 
989252845208 2 St i IS ER DC o “O TO 38 993928839 

кєз дсй م٥‎ ۵8 282858 م‎ ۹ δῶστοα о599:28390955090%49%9529 
3305931137400 820338755 5:9 72 6555 958559000 
مه هډ هم کح هه هډ ه لم م ٢ه هم به هڼڅ ېنځ په م‎ QZ که‎ ος هه‎ мыч = Кеч у. 


o مهم ےم‎ 
QAQNREDODÓSSOSESSIQSSS25020320535585929239235995 


12 31 42 62 84 34 71 
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© 


SSR 


4 ه۹۹۹ ۹۸ το‏ 
و ېم هوه همه هم هې هه o‏ 


рзд 


Ο 0 


86 
82 
82 
86 
1 86 52 73 14 53 82 ۱6١ 75 


ν 9 “5 
BRR OBS 
“44489 
5585348 


-- СҸ سس‎ 


5995958 


20 22 40 44 80 


9в888%Чс 
ras sss sa 


۳ه مه 8 


NAN 
Gg ga 


595548 


== ла ο em 


сум‏ با 00( م 
828982389 
5500544 


N — ο ο CY 
Ὃ مس زي‎ OO 


Ελ.‏ و وع 


ON ο ου 
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Chapter V 


Smarandache Concatenated Sequences 


V.1 Introduction 

Smarandache formulated a series of very artificially conceived sequences 
through concatenation. The sequences studied below are special cases of the 
Smarandache Concatenated S-sequence. 


Definition: Let G-(gi, £2, .... gx, .... ) be an ordered set of positive integers 
with a given property G. The corresponding concatenated S.G sequence is 
defined through 


S.G = {a,:a, = g,,a, =а, ,-10 هت‎ +9 к>] 


In table 1 the first 20 terms are listed for three cases, which we will deal with 
in some detail below. 


V.2 The Smarandache Odd Sequence 


The Smarandache Odd Sequence is generated by choosing 
G711,3,5,7,9,11,.....). Smarandache asks how many terms in this sequence 
are primes and as is often the case we have no answer. But for this and the 
other concatenated sequences we can take a look at a fairly large number of 
terms and see how frequently we find primes or potential primes. As in the 
case of prime-product sequence we will resort to Fermat’s little theorem to 
find all primes/pseudo-primes among the first 200 terms. If they are not too 
big we can then proceed to test if they are primes. For the Smarandache Odd 
Sequence there are only five cases which all were confirmed to be primes 
using the elliptic curve prime factorization program. In table 2 # is the term 
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number, L is the number of digits of N and N is a prime number member of 
the Smarandache Odd 
Tobie 1a. The first 20 terms of the Smarandache Odd Sequence 


























1 
13 
135 
1357 
13579 
1357911 
135791113 
13579111315 
1357911131517 
135791113151719 
13579111315171921 
1357911131517192123 
135791113151719212325 
13579111315171921232527 
1357911131517192123252729 
135791113151719212325272931 
13579111315171921232527293133 
1357911131517192123252729313335 
135791113151719212325272931333537 
13579111315171921232527293133353739 
1357911131517192123252729313335373941 


Table 2. Prime numbers in the Smarandache Odd Sequence 





# L N 

2 2 13 

10 15  135791113151719 

16 27 135791113151719212325272931 

34 63  1357911131517192123252729313335373941434547495153555 
75961636567 

49 93 1357911131517192123252729313335373941434547495153555 
75961636567697173757779818385878991939597 





76 


Term #201 is а 548 digit number. 
V.3 The Smarandache Even Sequence 


The Smarandache Even Sequence is generated by choosing G={2,4,6,8, 10, 
......}. Тһе question here is : How many terms аге nth powers of a positive 
integer? 


A term which is a nth power must be of the form 27.32 where a is an odd nth 
power. The first step is therefore to find the highest power of 2 which divides 
a given member of the sequence, i.c. to determine n and at the same time we 
will find a. We then have to test if a is a nth power. The Ubasic program 
below has been implemented for the first 200 terms of the sequence. No nth 
powers were found. 


Table 1b. The first 20 terms of the Smarandache Even Sequence 


2 
24 
246 
2468 
246810 
24681012 
2468101214 
246810121416 
24681012141618 
2468 101214161820 
246810121416182022 
2468 1012141618202224 
2468 101214161820222426 
246810121416182022242628 
2468 1012141618202224262830 
2468 101214161820222426283032 
2468 10121416182022242628303234 
24681012141618202224262830323436 
2468 101214161820222426283032343638 





246810121416182022242628303234363840 
24681012141618202224262830323436384042 
Ubasic program: {only the essential part of the program is listed) 


60 N=2 

70 for U%=4 to 400 step 2 

80 DR=int(log(U%) /log(10}}+1 ‘Determine length οἱ U% 
90 М-М”10А0%-/% ‘Add on U% 

100 A=N-E%=0 

110 repeat 

120 Al=A:A=A\2:inc EX ‘Determine EX (=n) 

130 until res<>0 

132 dec РЖ:А=АЇ ‘Determine A (=o) 

140 B-round(A^(T/EX)) 

150 if BAEX=A then print 1 ‘Check if a is anf power 
160 next 

170 end 


So there is not even a perfect square among the first 200 terms of the 
Smarandache Even Sequence. Are there terms in this sequence which are 2-p 
where p is a prime (or pseudo prime). With a small change in the program 
used for the Smarandache Odd Sequence we can easily find out. Strangely 
enough not a single term was found to be of the form 2-р. 


V.4 The Smarandache Prime Sequence 


The Smarandache Prime Sequence is generated by {2,3,5,7,11, ...}. Again 
we ask: - How many are primes? - and again we apply the method of finding 
the number of primes/pseudo primes among the first 200 terms. 


There are only 4 cases to consider: Terms #2 and #4 are primes, namely 23 
and 2357. The other two cases are: term #128 which is a 355 digit number 
and term #174 which is a 499 digit number. 


#128 

23571 1131 7192329313741 4347535961 6771 7379838997 101 1031071091 13127131 137139 
149153 1571631 67173179183 191 19319719921 |} 22322722923323924 | 2512572632692712 
7728128329330731 131331 7331 33734734935335 9367373379 383389 39 7 401 409419421 43 


1433439443449457461463467479487491499503509521523541547557563569571577587 
593599601607613617619631641643647653559661673677683691701709719 


Table Ic. The first 20 terms of the Smarandache Prime Sequence 






2 
23 
235 
2357 
235711 
23571113 
2357111317 
235711131719 
23571113171923 
23571 11317192329 
235711131719232931 
23571113171923293137 
2357111317192329313741 
235711131719232931374143 
23571113171923293137414347 
2357111317192329313741434753 
235711131719232931374143475359 
23571113171923293137414347535961 
2357111317192329313741434753596167 
235711131719232931374143475359616771 
23571113171923293137414347535961677173 






















#174 

23571 1 1317192329313741 4347535961677 17379838997 101 1031071091 13127131137139 
1491311571631671731 79181 1911931971992] 1223227229233239241 251 2572632692712 
7728128329 3307311313317331337347349353359367373379 383389397 401 40941942143 
143343944344945746 | 463467479487 49 1 499503509521 523541547557563569571577587 
993599601 60761361 7619631641 64364765365966 1 67367768369 1 701 7097 197277337397 
43751757761 76977378779780981 | 82182382782983985385 785986387 788 | 8858879079] 
1919929937941 947953967971 97798399 1 997 1009101310191021 10311033 


Are these two numbers prime numbers? 
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ПИЕ ποπ ομως 


Chapter VI 


On the Harmonic Series 


VL1 Comparisoa of a few sequences 

The harmonic series 14-545... is so important and so intensely 
studied that nothing more can be said about it - or can there? We may always 
pose a few questions. 


We know that the series is divergent. Furthermore the series Ils 
=1 


convergent for a>1 and divergent for as1. Since it is the “borderline” it may 
be interesting to examine the following questions: 

Question 1. What is the smallest number of terms m we need in order to 
make the sum greater than a given positive integer w? In other words the 
smallest m for which Στον 


Question 2. What is the smallest m for which زرو کنا‎ 
η бернн ουσ 
introduced in chapter П.2,3 we pose the corresponding questions: 


Question 3. What is the smallest m for which لر‎ w? 
1 


Question 4. What is the smallest m for which n iw? w? 
-1 


Table 1 shows the values of πι required in each case for w=1, 2, 3, ...12. In 
diagram 1 In(m) is plotted against w. سه‎ ο 


shown on the cover. We see that the series Στα is very “sensitive” to 
k-ik 


Changes in the value of a. ا‎ zu 
k=1 Š, 


a | : А 
Ж and > —— on the divergent side. 
k=1 S? (k) 5 


Table 1. Comparison of sequences. 


is on the convergent side of the graph for 


Series 1: Series 2: Series 3: Series 4: 
ы 1/k 1/s?(k) 1/S,*(k) 1/k}-° 
Values of m Values of m Values of m Values of m 


دم 
N‏ 
N‏ 
e‏ 
N‏ 
N‏ 


2 4 144 9 4 
3 11 462 54 12 
4 31 1045 243 35 
5 83 1944 729 102 
6 227 3200 2048 311 
7 616 4862 6561 966 
8 1674 6912 16384 3082 
9 4550 9477 37746 10115 
10 12367 12480 93312 34167 
11 33617 16008 209952 118971 
12 91380 20020 497664 427698 





Diagram 1. The behaviour of some sequences: тт] plotted against w. 
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VL2 Integers represented as sums of terms of the harmonic series 


Given a positive integer w. Can we represent w as a sum of a finite number of 
different terms chosen from the harmonic series? For w-1 we have the simple 


answer 1= eeu but what about w=2, w=3 etc. 


Method 1: k-perfect numbers. By definition a k-perfect number 15 a number 
n which satisfies the equation o(n)=k-n, where o(n) is the sum of divisors of 


(1 and n included). دادے‎ therefore gives the desired representation of the 


integer k. The result of a computer search for k-perfect numbers is given in 
table 2. A column is included for the number of divisors m since this is equal 
to the number of terms required for the representation. The numbers are grven 
in factor form as this is more interesting than the decimal representation. A 2- 
perfect number is in the literature simply referred to а perfect number. They 
have attracted a lot of interest - The problem of perfect numbers, a favorite 
with ancient Greeks, owes its origin to the number mysticism of the 
Pythagoreans (quote from Elementary Number Theory, Uspensky Heaslet). 


Table 2. k-perfect numbers п, m = the number of divisors. 








k = n 
2 4 2-3=6 
2 6 22.7-28 
2 10 231-496 
2 14 2°.127=8128 
2 26 212.8191=33550336 
2 34 216.131071-8589869056 
2 38 2*8.524287=137438691328 
22 62 2°°.2147483647=2305843008139952128 
3 16 23-3-5-120 
3 24 2°-3-7=672 
3 288 2°.5-7-19-37-73=459818240 
3 80 2°.3-11-31=523776 
3 224 213.3-11-43-127-<1476304896 
3 480 214.5.7-19-31-151-51001180160 


Table 2. continued. 


uu u-————t€—H————————————— U 2717 n rr 777777: 2:727 ———— 
ж n 


k 
4 216 22.32.5.72.13.19=2178540 
4 96 2?.32.5.7.13=32760 
4 96 2°3°-5-7=30240 
4 384 2733.52.17.31-45532800 
4 320 2?.37.5.11.31-23569920 
4 480 2?-3°-7-11-13-31=142990848 
4 1056 219.33.52.23.31.89-43861478400 
4 998: د کح ېږل‎ 32.52.19.31.683:2731.8 :8191 AA ee E 
5 1920 2/.35.5.7-117.17-19-14182439040 
5 2304 2'.3°.5.72.13.17:19=31998395520 
5 5280 219.3*.5.7.112.19-23-89 
5 6336 219.35.5.7?.13.19.23.89 
5 3456 211.33.54.12.13.19.31 
5 9216 211.35.63.73.133.19 
5 43008 2753.5 713517327337 
$ 709632 221.36.52.7.19.232.31.79-89-137-547-683-1093 
5 94208  2:5.3'-5-7.11-19-41-47-151-197-178481 
6 


тәеететететететтетететететтетатаеттетететеоттсттеттестетееее 


71658880 21731.52751]213219131:37.6173181 


For w=3 we һауе a representation by 16 terms: 

1 1 1 1 1 1 1 1 1 1 l 1 1 1 1 
З=1+—+—+—+—+—+—+—+—+—+—+—+—+—+-—+ 

2 3 4 5 6 8 10 12 15 20 24 30 40 60 120 
The smallest number of terms required for w=5 is 1920 and for 6 a staggering 
1,658,880 terms. For w-7 we have no representation at all. We would first 


have to find 7-perfect number. 


No great results. It's like having used a sledge hammer to kill a mosquito and 
missed the target. But k-perfect number are interesting in themselves. 


Method 2: Trivial Expansions. A term cin the harmonic series can always 
be replaced by two other different terms from the harmonic series: 


1 Е : . : 
~= = —— +————. Given an expansion of an integer w into terms 
κ k(k*]) К+1 k(k+b Р 
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of the harmonic series it is possible to construct infinitely many such trivial 
expansions of w. If however, in a given expansion of w, each term is replaced 
by two other terms as above and furthermore the replacement process is 
carried out until all terms in the new representation are different from one 
another and also from the terms in the original representation then a new 
representation has been obtained. Obviously such representations for w; and 
у, with no common terms can be added together to form a representation of 
W=W) FW. 


Example: 
ЕЕЕ 
2 3 
Expand: k > k+1 k(k*1) 
2 — 3 6 Not to be used (repeated) 
3 E 4 12 
6 — 7 42 
3 and 6 аге repeated and since their expansion gives 4, 12 and 7, 42 we have 
to continue the expansion: 
4 چ‎ 5 20 
12 — 13 156 
7 -» 8 56 


42 > 43 1806 
From this we can write the following representation of 2: 

1 و1‎ l1 i- 1 1 1 1 1 1 1 
2=1+—+—+—+—+—-+—+—+—-+—+—+—-+—— 

4 5 7 8 12 13 20 42 43 56 1% 6 
Adding our original representation of 1 we get the following representation 3: 

| 1 1 1 1 1 1 1 1 1 1 1 1 1 4 
3=1+—+—+—+—+—+—+—+—+—+—+— -+—-+— +— +—— 

2 3 4 5 6 7 8 12 13 20 42 43 56 1S 1806 
As in the case of 3-perfect numbers we have a representation of 3 with 16 
terms. The minimum number of terms for which a representation is possible 
is 11 as can be seen from table 1. Let's try to get closer to this. 


Method 3. Maximum density. For a given value of w table 1 gives the 
minimum number of terms from the harmonic series required to form a 
representation. For a given value of w=3 this number is m. We can therefore 
try to choose the first m-1 terms as close together as possible, i.e. we form 


where 5 is a small fraction which we һауе to expand іп harmonic terms. As 
before we apply the idea to w=3. We get 

1 1 1 1 1 1 1 1 1 1 1 1 

3=1+—+—+—+—+—+—+—+—+—+—+—+— 

2 3 4 5 6 7 8 9 10 15 230 57960 
Bravo! Only 13 terms. We had 16 with the other two methods. But no 
pleasure lasts forever. The nasty little fraction s will give no end to problem 
when we try large values for w. Although small s will have large numerators 
and denominators as we shall soon see. Maybe better get the sledge hammer 
out again and go in search for the 7-perfect number. 


VL3 Partial sums of the harmonic series as rational numbers. 


We consider consecutive partial sums S, = Ў = expressed as rational 
k=1 


numbers reduced to lowest terms. 


Question 1. Can there be a repetition of numerators? 
Assume 5,1 => and S, E where (a,b)=1 and (a,bj)-1. We then have 
1 


= Te which is equivalent to an(b+;), but (a,bb,)=1 and the latter 
1 


equation is therefore impossible and the assumption is wrong. We conclude: 


Two consecutive sums S,; and S, cannot have equal numerators when 
reduced to lowest terms. 


A necessary condition for fidig two sums S, = Е and S, = m , т>п+1 with 
n m 


equal numerators a,=a,, when in lowest terms is therefore that πι exists so that 
n<n <m and a, <a, |. Such numerators exist. Examples: 

ас-49 

815137 

ais=14274301 

а17=4214223 
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هوو وغو Ка‏ 


a99=360968703235711654233892612988250163157207 
81007 14466636279520351160221518043104131447711 
etc 


Now, let p be a prime then .هنېد‎ because S, is a monotonously increasing 
function of n and 


Conjecture: Let p; and pj: denote consecutive primes. The numerators in the 
interval p;<n<p,,,, after reducing to lowest terms, are all greater than a, —1, 


i.e. when a term of the type = is added the numerator takes a leap to a higher 


value from which, although it may decrease, it will never drop to a value 
below the value before the leap. 


Question 2. What is the most frequently occurring denominator in the 
sequence ος. ë + when reduced to lowest terms? 





Let 5,1-5 where (ab)-l. Then سوت و‎ s, Consider the 


greatest common divisor of b and n, d-(b.n). Then it is seen immediately that 
the fraction can be reduced by d. Put b-bid and n=n,d and consequently 
_ m +b, 
Sa = bnd 
If a further reduction is possible then it must be by a factor in d or by d itself 
because (п, ΡΞ! and (Οι ΑΕΙ. 


If n=p is a prime number then (b,n)-1 and the denominator os S, will have p 
as a factor. The denominator of S, is therefore unequal to all previous 
denominators and is at most the first in a series of equal denominators. Sums 
with equal denominators must therefore be contained between 


1 1 1 l 1 1 
5. =1+—+—+... +— and S, =1+=+=+... + 
b dec aan "23 Pj+ı 





where p; and pj; are consecutive primes. In table 3 the number of equal 
denominators is denoted q and the corresponding number of digits of the 
denominator L. 


Table 3. Equal denominators for n«1000. 





p Бат” а L 
89 96 8 39 
1:39 148 10 62 
3117 330 14 140 
891 906 16 389 





I hope this book has been more meaningful than reading the number below 
which is the 389-digit denominator which occurs 16 times in consecutive 
partial sums of the harmonic series for n<1000. 


4179982336319706196962068134998944512773598562712 
1827499073959690180767635682039557104883261454628 
9717289725066739371544332075631263681868 936793360 
4957159732797690125645682733080622438318699474351 
6944805893417522011858358480060385932366764886574 
4940462245050719348300672374671465097882430105554 
5638386382683543980298959618721940328647837326571 
6039813081070459435310984219943740583597120000 


